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Abstract 



This paper is a continuation of an earlier work (arXiv:1210.7928) which computed 
analytically the tree-level annihilation rates of a collection of non-relativistic neu- 
tralino and chargino two-particle states in the general MSSM. Here we extend 
the results by providing the next-to-next-to-leading order corrections to the rates 
in the non-relativistic expansion in momenta and mass differences, which include 
leading P-wave effects, in analytic form. The results are a necessary input for the 
calculation of the Sommerfeld-enhanced dark matter annihilation rates including 
short-distance corrections at next-to-next-to-leading order in the non-relativistic 
expansion in the general MSSM with neutralino LSP. 



1 Introduction 



The increasing precision on the experimental determination of the dark matter (DM) 
density, which is expected to be further improved by the data from the PLANCK satellite, 
has brought a renewed interest in the impact of radiative corrections to the annihilation 
cross section of dark matter candidates of particle nature. Particles with weak interaction 
strength and masses around the TeV scale that dropped out of thermal equilibrium in the 
Early Universe yield the correct order of magnitude for the relic density. An example of 
the latter is provided by the lightest neutralino of the minimal supersymmetric extension 
of the standard model (MSSM), perhaps the most promising candidate for such weakly 
interacting dark matter. The potential to set stringent constraints on the parameter 
space of the MSSM using the high precision measurements of the dark matter density 
crucially depends on having an accurate calculation of the neutralino relic abundance. 

A necessary input for this calculation is the annihilation cross section of the lightest 
neutralino, and of all possible co-annihilation processes. While programs exist that 
provide the tree-level results numerically [HE], at the one-loop level such calculations 
are not available for a generic MSSM model, although they have been performed for 
some scenarios [2H7], or under certain approximations [HHTU]. 

There is however a certain class of radiative corrections where higher-order loop 
diagrams contributing to the DM annihilation amplitude are not necessarily suppressed. 
At the temperatures where freeze-out of the relic particle abundance takes place, the dark 
matter particles are non-relativistic, with typical velocities of order v ~ 0.2 c. Quantum 
loop corrections due to the exchange of light particles between the non-relativistic DM 
particles before annihilation can become more and more important in situations where 
the force coupling strength is larger than the DM velocity and the mass of the force 
carrier is much lighter than the DM mass, eventually requiring a resummation of the 
terms in the perturbative expansion to all loop orders. This phenomenon has been 
termed as "Sommerfeld effect", and can lead to a significant enhancement of the DM 
annihilation rates, also of relevance for the calculation of primary decay spectra in the 
present Universe. In the MSSM, Sommerfeld corrections may constitute the dominant 
radiative correction when the lightest neutralino is much heavier than the electroweak 
gauge bosons and the Yukawa potential generated by their exchange becomes long-range. 
For such heavy neutralinos mixing effects are suppressed by 0(Mz /^lsp) and thus mass 
degeneracies arise in the neutralino-chargino sector, making necessary to account for co- 
annihilations processes in the relic density calculation. Two prominent examples of this 
scenario are the MSSM wino- and Higgsino-limit, for which the impact of the Sommerfeld 
effect has been extensively studied [TTHT5] 

The tree-level DM annihilation cross section can be expanded in the relative velocity 
v TC \ of the two annihilating particles (v ie \ = \v\ —1)2]), 

ffann ^rel = d + b + 0(v* el ) , (1) 

for non-relativistic v Te \. In a previous work [16] (referred to as paper I in the following) 
the calculation of the leading-order coefficient a in (JTJ was presented in analytic form for 
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the general MSSM with neutralino LSP, including results for all co-annihilation processes 
with nearly mass-degenerate neutralinos and charginos. In this paper we complete the 
calculation by providing analytic results for the subleading term in this expansion, the 
coefficient b. In the non-relativist ic effective theory (EFT) framework devised in paper 
I, the coefficients a, b account for the short-distance part of the neutralino and chargino 
pair-annihilation processes, and are written as a combination of the absorptive parts 
of Wilson coefficients of local four-fermion operators. The absorptive part of the EFT 
matrix element of these four-fermion operators then gives the full neutralino and chargino 
pair-annihilation rates, including the long-range Sommerfeld effect. 

While the leading-order term a in (CQ) receives contributions only from S-wave anni- 
hilations, the subleading term b encodes both S- and P-wave annihilation contributions, 
which we provide separately. This is required for a correct implementation of the Som- 
merfeld correction factors, that depend on the spin and partial- wave configuration of the 
annihilating state. In recent literature which addresses the Sommerfeld effect includ- 
ing the 0(v^ el ) terms in the annihilation cross section, only the P-wave contributions 
to the coefficient b have been computed approximately using numerical routines at the 
amplitude level [ID] , or it has been assumed that b is entirely P-wave [17] . 

The results presented in I (and complemented in this work) also extend those from 
previous approaches in another relevant aspect. The knowledge of the tree- level annihila- 
tion cross section (Jl]) for all nearly mass-degenerate neutralino and chargino two-particle 
states is not sufficient for the calculation of the Sommerfeld corrected (co-) annihilation 
rates. As described in paper I, a contribution to the full annihilation rate of an incoming 
XiXj state is given by the imaginary part of the amplitude for the process 

XiXj -»■ • • • Xe 1 Xe 2 -> X A X B -»■ XeiXe 3 ~* ■ ■ ■ ~> XiXj , (2) 

where the transitions among the XX states are mediated by long-range potential interac- 
tions and the short-distance annihilation into SM and light Higgs particles (XaXb) in- 
volves the two-particle states XeiXe 2 an d Xe 4 Xe 3 - The case when XeiXe 2 &nd Xe4Xe 3 dif- 
ferent states corresponds to an off-diagonal short- distance annihilation rate. In our EFT 
approach the diagonal and the off-diagonal (tree-level) short- distance rates are encoded in 
the absorptive part of local four-fermion operators' Wilson coefficients, that are obtained 
from matching the EFT tree-level matrix elements of the four-fermion operators to the 
absorptive part of the hard (1-loop) MSSM amplitudes for the Xe 1 Xe 2 ~^ XaX b — > Xe 4 Xe 3 
scattering reactions. The off-diagonal terms have not been taken into account in the 
Sommerfeld-enhanced neutralino relic abundance calculations aside from the wino- and 
Higgsino-limits [TTHT31 IT5]. and their implementation using the numerical packages that 
provide the tree- level annihilation rates has not yet been attempted. In contrast, the 
analytic results presented in this work allow for a systematic treatment of all diagonal 
and off-diagonal short- distance annihilation rates at next-to-next-to-leading order in the 
non-relativistic expansion in Sommerfeld-enhanced annihilation reactions. 

The contents of this paper are the following: In Sec. |2]we briefly review the effective 
Lagrangian framework introduced in paper I and recollect the essential notation. We 
then introduce the dimension-8 four-fermion operators that encode the next-to-next-to- 
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leading non-relativistic corrections to the short-distance annihilation of neutralino and 
chargino particle pairs. As for the case of the (leading-order) dimension-6 operators dis- 
cussed in paper I, the analytic results for the absorptive part of the Wilson coefficients 
can be obtained as the product of kinematic and coupling factors. In the Appendix [A] 
we recall the master formula to write down the Wilson coefficients and the rules for its 
implementation. While we rely on paper I for the extraction of coupling factors also 
for the Wilson coefficients presented in this work, explicit expressions for the P-wave 
kinematic factors are given in the Appendix |A] The expressions for the (rather lengthy) 
next-to-next-to- leading S'-wave Wilson coefficients are collected in a Mathematica pack- 
age attached to this paper [18], which also includes the P-wave kinematic factors and 
those from the leading-order operators that were written explicitly in the appendix of 
paper I. Appendix |B1 explains the notation used in this electronic supplement. In Sec. |3] 
we generalise the formula for the tree-level annihilation rates to the case of off-diagonal 
annihilation processes, which is needed in order to analyse the size of next-to-next-to- 
leading corrections in an off-diagonal transition used as a case example in Sec. HI Apart 
from the latter, Sec. H] also discusses three selected (diagonal) processes where the role of 
the next-to-next-to-leading corrections from our analytic calculation of the annihilation 
cross section is markedly different. Through these examples we illustrate the importance 
of separating the different partial-wave contributions to the short-distance annihilation 
for the computation of the Sommerfeld-corrected cross sections. Our results for the 
diagonal annihilation rates for these examples are checked against the corresponding 
unexpanded cross sections computed with a numerical code. Finally we summarise our 
findings in Sec. [5j 



2 Basis of the dimension 8 operators in ££ a nn 

The non-relativistic MSSM (NRMSSM) effective theory set-up of paper I is built out 
of no < 4 nearly on-shell non-relativistic neutralino (x°, i = l,...,no) and n + < 2 
nearly on-shell non-relativistic chargino (Xj,j — ^-i---i n +) modes whose masses are 
nearly degenerate with the mass tt^lsp of the lightest neutralino Xi- As pair-annihilation 
reactions of these neutralino and chargino species into (not non-relativistic) SM and 
light Higgs-particle final states take place at distances of the order C(l/mLSp), much 
smaller than the characteristic range of potential interactions between the incoming 
XX two-particle states, we can incorporate the actual annihilation rates in the effective 
theory through the absorptive part of Wilson coefficients of local four-fermion operators 
(5£ ann ), in analogy to the treatment of quarkonium annihilation in NRQCD [TJJ]. In 
contrast to the QQ case, here the long-range potential interactions can lead to transitions 
among different two-particle states before the short-distance annihilation takes place: the 
initially incoming XiXj particle pair can scatter to any accessible (nearly on-shell) Xe a Xe b 
particle pair prior to the annihilation into SM and light Higgs two-particle final states 
XaX b . Therefore we have to account for absorptive parts of generic XeiXe 2 — * XaX b — > 
Xe 4 Xe 3 amplitudes, where XeiXe 2 an d Xe 4 Xe 3 can be different states (see Fig. 1 in paper 
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I)- 

At 0(a 2 ), where a = g 2 /An with g the SU(2) L gauge coupling in the MSSM, the 
absorptive part of the four-fermion operators' Wilson coefficients are obtained by match- 
ing the absorptive part of XeiXe 2 ~ * Xe 4 Xe 3 1-loop scattering amplitudes evaluated in the 
MSSM with the tree-level matrix element of four-fermion operators contained in 5£ a nn 
in the effective theory. At the 1-loop level, the contribution to the absorptive part from 
every individual final state X^Xb is free from infrared divergences, and can be given 
separately. At higher-orders the absorptive part of the Wilson coefficients refers to the 
inclusive case, i.e. summed over all accessible final states. 

The leading-order contributions to 5C ann are given by dimension-6 four-fermion oper- 
ators, encoding leading-order S"-wave scattering reactions. The corresponding operators 
have been given in paper I. At next-to-next-to leading order in the non-relativistic expan- 
sion in momenta and mass differences, dimension-8 four-fermion operators contribute^ 
Here we adopt the same notation used for in I in order to write the dimension-8 

four-fermion operators in SC^^ as 



5C 



d=8 
ann 



E l fXX^XX (lp \ QXX^XX lip \ 

AM 2 { eie2 H e4e 3} ^ 1 ' {e4e3}{e 2 ei} I l ) 

xx-^xx 



_j_ 1 fXX^XX (3p \ QXX^XX (3p \ 

Z 4 ' 4 A l\/f 2 { e l e 2}{e 4 e 3 } V J ) {e4e 3 }{e 2 ei} V J ) 



A^M 2 { e i e 2}{ e 4e3> 
XX^XX .7=0,1,2 



i_ . „XX-^XX (2s+l o \ -pXX^XX f2s+l rr \ 

+ Z^ Z^ A M2 9 {^}{e A e 3 } \ 1 {e 4 e 3 }{e 2e i} I 



AM 2 

xx->xx 8=0,1 



I V V - /iXX-i-XX (2s+ln\ ClXX^rXX (2s+\ q \ fo\ 

Z^ Z^ Z^ i W2 'Meie 2 }{e 4 e 3 } I °s) {e 4 e 3 }{e 2 e^ I °s) " W 



AM 2 

XX-^XX s=0,l 1=1,2 

The label XX ~~ >* XX stands for all (off-) diagonal non-relativistic scattering processes 
among neutralino and chargino two-particle states. Neutral reactions involve x°X° an d 
X~X + states, while singly-charged and doubly-charged processes include x°X ± an d X ± X ± 

states, respectively. The /g^JJU,}' ^H^} and ^K}{e 4 e 3 } in © denote the 
Wilson coefficients of the corresponding four-fermion operators 0{ eie3 y{ e2ei y, V{ eie3 }{e 2 e 1 } 
and Qi{e 4 e 3 }{e2ei}, whose explicit form for the case of XeiXe 2 — Xe 4 Xe 3 scattering reactions 



is given in Tab. [To The labels in (j3J) range over ej = 1, . . . , n (e» = 1, . . . , if 
the respective field Xe 4 in the XeiXe 2 — >* Xe 4 Xe 3 reaction refers to neutralino- (chargino-) 



1 Let us remark that we do not consider next-to- leading order contributions to <5£ ann , corresponding 
to dimension-7 four-fermion operators, as these encode 1 Sq — 3 Pq, 3 Si — 1 P\ and 3 Si — 3 P\ transitions 
which will require the addition of u rc i-supprcssed potential interactions in the long-range part of the 
annihilation (we consider only 0(v^ cl ) effects from the short-distance annihilation, and not those arising 
from sub-leading non-Coulomb (non- Yukawa) potentials, in consistency with paper I). 

2 In order to ensure the f7(l) em gauge invariance of the NRMSSM, all derivatives d in dimension-8 
four-fermion operators O and V that act on chargino fields (j\%,C,i) have to be replaced by the corre- 
sponding covariant derivative D = d + i eA, where A denotes the spatial components of the photon 
field A^. 
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QXX->-XX(± p.jj 

QXX^XX(3 p i > 
£)XX^XX(3p 2 ) 



Se. 



C4 



Si 



C J 2 



e e 1 (- 



2 ^e 4 



x«r)£ 



(3 



t'2 



^e 4 ^ 2 " 



Si 



(-2 



-I ^)e ei 



-J V(v^) 



e el 



pxx^xx(i,5 ) 

pXX^XX(3 5l ) 



Se 4 Se 3 



■5 a) £ ei +£, 



t^\ 2 



Se 4 " ^e 3 



*e 2 



Qf c_ ' xx ( 3 5 , i) 
g xx-xx ( icy 



Table 1: Explicit form of the P-wave (O) and next-to-next-to leading order S"-wave (V, 
Qi) four-fermion operators contributing to XeiXe 2 — Xe 4 Xe 3 scattering reactions. Each 
index e, can take the values e, = 1, . . . , no- The P- and next-to-next-to-leading order 
S'-wave four-fermion operators for the remaining neutral, charged and double-charged 
Xe 1 Xe 2 Xe 4 Xe 3 processes are obtained by replacing the field operators £ ei , i = 1, ... ,4 
above by those of the respective particle species involved. The quantity d is a 3-vector 
whose components are d l = d/dxi. The action of d on the two field operators at its 
left and right is defined as d £ eo = $ (d£ ea ) — (d^ Y £ Ca . The symmetric traceless 
components of a tensor T ij are denoted by = (T ij + T ji )/2 - T kk 5 ij /3. Finally, the 
mass scale M is defined in (j4j) and the mass differences 5m, 5m are given in (JSJ). 



species. The factor 1/4 in front of the operators in ([3]) is a convenient normalisation of 
transitions matrix elements in the effective theory. In addition, a normalisation factor of 
1/M 2 has been factored out in ([3]), such that the next-to- next-to-leading order Wilson 
coefficients have the same mass-dimension (—2) as the leading-order ones presented in 
I. The mass scale M is equal to half the sum of the masses of the x ei particles involved 
in the reaction x ei Xe 2 ~> Xe 4 Xe 3 , »-e. 

I 4 

M = 2 S m ^ ' ^ 

i=l 

such that M itself constitutes a process specific quantity. The quantum numbers 2s+1 Lj 
of the operators in SCf^n correspond to the angular-momentum configuration of the an- 
nihilating two-particle state. Note that the operators Qi ( 2s+1 S s ) have the same structure 
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as the dimension-6 operators O ( 2s+1 S s ) defined in paper I, but are proportional to the 
mass differences 

= m e4 -m ei = m e3 - m e2 

2 2 

computed from the masses m e . in the reaction XeiXe 2 ~~ >" Xe 4 Xe z - The mass differences 
§5§ have to be considered as 0(v 2 el ) effects in the expansion of the amplitudes according 
to the discussion given in Sec. 2.4 in I. Since 5m = 5m = for diagonal annihilation 
reactions Xe 1 Xe 2 — > Xe 1 Xe 2 (where the absorptive parts of the respective amplitudes are 
related to the corresponding annihilation cross section), the Qi ( 2s+1 S s ) are only relevant 
for the computation of the off-diagonal rates. 

We note that dimension-8 operators P( 3 5'i, 3 Di), which describe 3 Si — > 3 D± transi- 
tions, have not been included in 5C^. In the calculation of the tree- level annihilation 
cross section in the centre-of-mass frame, contributions from these operators vanish, 
while for the Sommerfeld enhanced annihilation cross section they will require to con- 
sider a ^-suppressed potential interaction in the long-range part of the annihilation in 
order to compensate for the change in orbital angular momentum in the short- distance 
part, thus yielding a contribution to the cross section of 0(vf el ). 

As already discussed in I, we construct £ ann in such a way that it contains all redun- 
dant operators, which arise through interchanging the single-particle field-operators at 
the first and second (third and fourth) position given a specific four-fermion operator, 
such that several operators describe one specific scattering reaction with a Xei and Xe 2 
(Xe 4 and Xe 3 ) particle in the initial (final) state. The respective Wilson coefficients re- 
flect the redundancy in symmetry relations under the exchange of the respective particle 
labels. Generalising from the leading-order 5"-wave relations given in Eq. (8) in I, the 
relations read in case of Wilson coefficients associated with operators O and V in ([3]) 

jX e2 Xe x ^Xe 4 Xe 3 (2s + 1t \ _ /_1\S + L jJCej Xe 2 ~^Xe 4 Xe 3 (2s + l T \ 

{e 2 ei}{e 4 e 3 } I J ) \ ) {eie 2 }{e 4 e 3 } I J) ' 

k X B1 X^Xe 3 X B4 (2s+X Lj S = r^s+L k X ei Xe 2 ^X e4 Xe 3 (2.+X L \ (6) 
{eie 2 }{e3e 4 j \ J l \ 1 \e\e 2 \\e 4 ez) \ J J ' V / 

where k — f,g for P- and next-to-next-to-leading order S'-wave coefficients, respectively. 
Finally note, that the hermiticity property of the non-relativistic Lagrangian leads to 
the relation 

UXX^XX (2s+1t \ _ \i,XX->XX (2s+1t \~|* f 7 \ 

K {eie 2 }{e 4 e 3 } { ^J)~ [ K {e 4 e 3 }{ ei e 2 } \ V ) 

for Wilson-coefficients k = f,g associated with the P- and next-to-next-to-leading order 
S'-wave operators O and V, in analogy to the respective relation given in Eq. (13) in I. 
Similar relations as flHJ Cj) above apply for the Wilson coefficients hi, where however an 
additional exchange of the particles in the definition of the mass differences 5m, 5m in 
front of the corresponding operators Qi has to be taken into account. 
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3 (Off-) diagonal annihilation rates r Xe i Xe 2^ XAXB ^ Xe * Xe 3 

In order to assess the importance of the non-relativistic corrections computed in this 
work, we shall compare in Sec. H] the EFT and full results for the tree-level annihilation 
rates for some selected processes. To extend this analysis to the case of off-diagonal 
annihilation rates, we generalise the definition of the rates in the following way. We 
define the centre-of-mass frame tree-level annihilation rate T associated with the (off-) 
diagonal Xe 1 Xe 2 ~ * XaXb —> Xe 4 Xe 3 scattering reaction as the product of the Xe 1 Xe 2 ~^ 
XaXb tree-level annihilation amplitude with the complex conjugate of the tree-level 
amplitude for the Xe^Xes ~* XaXb annihilation reaction, integrated over the final XaXb 
particles' phase spacqj and averaged over the spin states of the respective incoming 
particles Xe t , i = 1, • • • , 4. In the latter spin-average it is assumed that the XeiXe 2 an d 
Xe 4 Xe 3 pair reside in the same spin state@ In terms of the Wilson coefficients of the four- 
fermion operators in 5£ ann , the expansion of the annihilation rate V in the non-relativistic 
momenta and in the mass differences 8m, 5m is then given by 

+ (f( l Pi) + ~ f( 3 Po) + f( 3 Pi) + 1 h 3 P2)) 

+ 4^(^ o)+3 ^ (35i) ) 

+ 0((p 2 + p' 2 ) 2 , (p-p 1 ) 2 , p {,)2 5m, p^ 2 5m, 5m5m) , 

where p and p' correspond to the momenta of the \ ex and Xe 4 particle, respectively, in 
the centre-of-mass frame of the reaction. To shorten the notation we have suppressed 
in the label "XeiXe 2 ~^ XaX b — > Xe 4 Xe 3 " on the Wilson coefficients f,g and hi. As 
we study annihilation rates of non-relativistic Xe a Xe b particle pairs, the mass differences 
5m and 5m have to be (at most) of the order of the Xe a Xe b non-relativistic kinetic 
energy, as argued in Sec. 2.4 of paper I. Note that the non-relativistic expansion (jSJ) 
incorporates this convention and assumes that 5m, 5m ~ 0(p 2 /M). In case of diagonal 
XeiXe 2 XaX b —> XeiXe 2 scattering reactions, the definition of the corresponding 

3 The product of tree-level annihilation amplitudes has to be multiplied with an additional symmetry 
factor of 1/2 if the final state particles are identical, Xa = Xb- 

4 In the calculation of Sommcrfeld enhanced XiXj ~ * XaXb pair-annihilation rates through the imag- 
inary part of the XiXj —>•••—> XeiXe 2 ~~ ► XaXb — > Xe 4 Xe 3 —>•••—> XiXj forward scattering reaction, 
the assumption that the incoming and outgoing particle pairs in the XeiXe 2 ~~ ► XaXb — > Xe 4 Xe 3 short- 
distance annihilation part have the same spin state implies that just leading-order potential interactions 
in the XiXj —>•••—> XeiXe 2 an d XiXj Xe 4 Xe 3 scattering reactions are considered, since the 

long-range potentials are spin-diagonal only at leading order and hence pass the spin-configuration of 
the incoming XiXj P air to tnc X ei Xe 2 and Xe 4 Xe 3 pairs. 
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annihilation rate T obviously coincides with the definition of the spin-averaged centre-of- 
mass frame tree- level XeiXe 2 ~* X^Xb annihilation cross section times relative velocity, 
a xe 1 xe 2 -^x A x B v ^ an j eX p ans i on j n (Jgj^ with p' = p, reduces to the non-relativistic 
expansion of a Xe ^ Xe ^ XAXB v TC \ as given in Eq. (27) of paper I (see also Eqs. (19T TT21 below). 



4 Results 

In paper I we have presented several examples for the numeric comparison of the non- 
relativistic approximation to the tree-level centre-of-mass frame annihilation cross-sec- 
tion <7 XeiXe2_>XAXs times relative velocity v ie \ = \v ei — v e2 \, 

aXeiXe2 ^X A X B Urei = a + b ^ + ( v 4j ^ (g) 

with the corresponding unexpanded result obtained with MadGraph [20]. The coeffi- 
cient a in the expansion ([9]) is expressed in terms of the leading order S'-wave Wilson 
coefficients as 

a =fCS ) + 3f( 3 S 1 ) , (10) 
and the coefficient b can be written as the sum b = bp + b$, where 

,2 



b 



S 



. fi'Pi) + 3 f(fPo) + fi'Pi) + 3 /(^ 2 )J , (11) 
(g( l S )+3g( 3 S 1 )) , (12) 



M 2 



and 

/W 2 = me f £2 (13) 
m ei + m e2 

is the reduced mass of the XeiXe 2 two-particle state. We have again suppressed in ( TTOl - 
IT21 the labels on the Wilson coefficients / and g that indicate the specific XeiXe 2 — >* 
XaX b — > XeiXe 2 reaction under consideration to simplify the notation. The prefactor 
(/i eie2 /M) 2 in front of the Wilson coefficients in f[TTHT2T) is needed to translate the cross 
section's expansion in p 2 , Eq. (jSJ) with p' — p, to the t> 2 el expansion used in (jHJ). 

From the comparison of the non-relativistic approximation to a XeiXe2 ^ XAXB v TC \ with 
the full result from MadGraph, it was shown in paper I that the non-relativistic ap- 
proximation reproduces the behaviour of the exact tree-level cross section times rela- 
tive velocity within a percent level deviation up to v rc \/c ~ 0.60 The numeric extrac- 
tion of the coefficients a and b from MadGraph data by means of a parabola fit to 



5 In case of processes with vanishing 5-wave contributions the agreement between MadGraph and 
the non-relativistic approximation is a bit worse, but still with an accuracy at the level of ~ 6% for 
v rcl /c ~ 0.4. 
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a xeiXe 2 -±x A x B Vrel j n non-relativistic regime thus provides a useful numeric check 
for the sum of leading-order S- wave Wilson coefficients f fTUj) . as well as for the sum of 
next-to-next-to- leading order S-wave and P-wave Wilson coefficients in f fTT|T2|) . How- 
ever, a splitting of the numerically extracted coefficients a and b into their constituting 
partial-wave contributions is not straightforward from publicly available numeric codes, 
as this requires manipulations at the amplitude level. The separate knowledge of the 
different 2s+1 Lj partial wave contributions to the tree-level (co-)annihilation rates is es- 
sential for a precise determination of Sommerfeld enhanced neutralino (co-)annihilation 
cross sections, because the Sommerfeld enhancements depend both on the spin- and or- 
bital angular momentum quantum numbers of the annihilating particle pair. Therefore 
a consistent treatment of the Sommerfeld enhancement including P-wave effects requires 
the separate knowledge of all relevant (off-) diagonal tree-level 1 Sq and 3 Si partial-wave 
annihilation rates both at leading and next-to-next-to-leading order, as well as the indi- 
vidual (off-) diagonal tree-level 1 Pi and 3 Pj partial-wave annihilation rates. In the latter 
case, the knowledge of the (spin- weighted) sum over the three different 3 Pq, 3 Pi and 3 P2 
partial-wave Wilson coefficients, 

f( 3 Pj) = IfCPo) + f( 3 Pi) + jj/( 3 P0 , (14) 

is sufficient, as long as only leading-order non-relativistic potential interactions between 
the neutralino and chargino states are taken into account in the full annihilation am- 
plitudes. This is because the leading-order potential interactions depend on the spin 
(S = 0, 1) of the Xe a Xe b particle pairs taking part in the XiXj XeiXe 2 ~~ ^ 

XaXb —> Xe 4 Xe 3 ... XiXj scattering process, but do not discriminate among the 
three spin-1 P-wave states 3 Pj with different total angular momentum J = 0, 1, 2. 

Recently, Sommerfeld corrections including P-wave effects have been subject of study 
at f-loop [10] and with full resummation [IT]. In these studies, the next-to- next-to- 
leading order contributions in the expansion of the relevant (co-) annihilation rates were 
assumed to be given only by P-waves. While such reasoning is justified when the leading- 
order S'-wave contributions to the annihilation rates are strongly suppressed with respect 
to the next-to-next-to- leading order coefficients in ([HI ED , it does not hold for the general 
case. In particular, P- and next-to-next-to-leading order S'-wave terms can come with 
differing signs, such that a partial compensation of different next-to-next-to-leading order 
contributions to the annihilation rates may occur. 

In order to illustrate the different behaviour of the P- and next-to-next-to-leading 
order S-wave contributions to the tree-level annihilation cross sections, we show in 
Figs. [1] and [2] results for the tree-level annihilation cross section times relative veloc- 
ity, a Xe ^ Xe ^ A B v re \, for three different processes. The plots refer to the same SUSY 
spectrum that was used in paper I, which contains a wino-like neutralino LSP with mass 
m x o = 2748.92 GeV, and an almost mass-degenerate wino-like chargino partner with 
m + = 2749.13 GeV. The next-to-lightest chargino state has a mass m + = 3073.31 GeV. 
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Figure 1: Numeric comparison of the non-relativistic approximation (solid lines) to the 
tree-level annihilation cross section times relative velocity, crv Te i, for the XiXi ~~ W + W~ 
process (left) and for the S-wave dominated reaction XiXi ~^ u d (right) with the corre- 
sponding unexpanded annihilation cross sections produced with MadGraph. Numeric 
errors on the MadGraph data are given by a v m \/ yfN, where N = 10 5 gives the num- 
ber of events used in the MadGraph calculation of each cross section value. The 
dash-dotted red (dashed black) curves represent the constant leading-order term in the 
non-relativistic expansion of the cross section plus the P-wave (next-to-next-to-leading 
order S'-wave) contribution, a + b P v^ cl (a + b s v^ el ). 

4.1 Example 1: XiXi — ► W + W 

The plot on the left hand side in Fig. [TJshows the XiXi ~^ W + W~ tree-level annihilation 
rate, a relevant co-annihilation rate in the neutralino LSP relic density computation. 
The solid blue line corresponds to the non-relativistic approximation to the tree-level 
annihilation cross section, a Xl Xl ^ w+w <y rel) and the points correspond to the full tree- 
level result obtained with MadGraph. The deviation between our approximation and 
the MadGraph data is at one percent level for v TC \/ c ~ 0.6 and in the permille regime for 
smaller relative velocities. Further, the composition of the non-relativistic approximation 
to cr Xl Xl ~* w+w v re \ out of P- and next-to-next-to-leading order S'-wave contributions 
can be read off from Fig. [I] the dash-dotted red line represents the contribution a+bp v^ el 
to OH]), while the dashed black line is a + bs v^ el . While both bp and bs are roughly of 
the same order of magnitude, the summed P-wave contributions enter with a positive 
sign (bp = 1.86 • 10~ 27 cms), whereas the summed next-to-next-to-leading order S-wave 
contributions come with a negative weight, bs = —0.88 • 10 _27 cms. It is worth to 
note, that the sum of next-to-next-to-leading order corrections in the XiXi ~^ W + W~ 
tree-level cross section times relative velocity gives a ~ 6% correction to the leading- 
order approximation for v TC \/c ~ 0.4. For this relative velocity, the corrections to the 
leading-order approximation from P- waves only amount to ~ 11%, while those from 
next-to-next-to-leading order S-wave contributions amount to ~ —5%. Hence, in the 
light of the expected future experimental precision on the measured dark matter density, 
it is crucial to take these corrections into account. Further, as generically the Sommerfeld 
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enhancements for each of the contributing partial waves are different, it will be needed to 
investigate the Sommerfeld enhanced annihilation cross section including P- and next-to- 
next-to-leading order S*-wave enhancements separately. This study within our formalism 
is postponed to [21] . 

The fact that the P-wave terms in the example of Fig. [1] contribute with positive 
sign is generic: the sum of all 2s+1 Pj partial- wave contributions to any XeiXe 2 — ► XaXb 
annihilation cross section has to be positive, as it results from the absolute square of 
the coefficient of the 0(p) terms in the expansion of the annihilation amplitude. More- 
over, the separate 2s+l Pj partial-wave contributions must also be positive, since different 
2s+1 Pj-wave amplitudes do not interfere in the absolute square of the annihilation am- 
plitude due to total angular-momentum conservation and the additional conservation of 
spin in the non-relativistic regime. The next-to-next-to-leading order S-wave contribu- 
tions to the Xe 1 Xe 2 XaXb annihilation cross section, however, result from the product 
of leading-order and next-to-next-to-leading order S- wave contributions in the expansion 
of the XeiXe 2 ~> XaX b amplitude. There is a priory no reason why this product should 
be positive, and hence negative next-to-next-to-leading order S-wave contributions to 
the cross section can occur, as can be explicitly seen in the examples presented in this 
section. 

4.2 Example 2: XiXi" ~~ ^ u d 

The right plot in Fig. [1] shows results for the S*-wave dominated tree-level XiXi" ~~ u< ^ 
annihilation process, also of importance in the neutralino relic abundance calculation 
including co-annihilations. The dashed black line, representing the a+bs w 2 el contribution 
to the non-relativistic expansion of the annihilation rate with b$ = —0.78 • 10~ 27 cms, 
basically coincides with the solid blue line, which corresponds to the complete non- 
relativistic approximation fl9]). Data produced with MadGraph for the XiXi ~~ >" u d tree- 
level annihilation rate are shown in addition, illustrating once again the nice agreement 
of the non-relativistic approximation with the unexpanded tree-level cross section results 
for relative velocities up to v rc \/c ~ 0.6. It is worth to understand the suppression of 
P-waves with respect to the next-to-next-to-leading order S'-wave contributions in the 
XiXi ~~ >* u d process as well as the composition of the coefficient bs out of its 1 S and 3 Si 
partial-wave contributions: First note, that in the case of vanishing final state masses, 
m u = rrtd = 0, the contributions to both a and bs can be attributed solely to 3 S\ 
partial waves. The absence (or more generally the suppression in m q /M, q = u,d) of 1 5'o 
partial- wave contributions both in the leading-order coefficient a and in bs is a helicity 
suppression effect. The helicity suppression argument applies to all 2s+1 Lj partial- wave 
reactions with J = 0, as the final state of a massless (left-handed) quark and a massless 
(right-handed) anti-quark in its centre-of-mass system cannot build a total angular- 
momentum state J = 0. Hence both 1 S'o as well as 3 Po partial- wave contributions are 
helicity suppressed. 

The suppression of P±, 3 Pi and 3 P2 partial- wave contributions that proceed through 
single s-channel gauge-boson or Higgs exchange is related to either factors of A m = 
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Figure 2: Left plot: Numeric comparison of the non-relativistic approximation (solid blue 
curve) to the tree-level annihilation cross section times relative velocity, crv Te \, for the 
P-wave dominated X2X2 — ^ h°h° reaction to data for the corresponding unexpanded an- 
nihilation cross section produced with MadGraph. Numeric errors on the MadGraph 
data are taken to be a v TC \/y/N, where N = 10 5 gives the number of events used in the 
MadGraph calculation of each cross section value. The dash-dotted red and dashed 
black lines represent the constant leading-order term plus the P-wave or the next-to- 
next-to-leading order S-wave contribution, a + bpv^ cl or a + bsv^ el , respectively. Note 
that the a + bpv^ el contribution and the non-relativistic approximation coincide, as there 
are no S-wave contributions in this particular annihilation reaction. Right plot: Off- 
diagonal annihilation rate T for the reaction xfXi ~^ W + W~ — > X2X2 ■ The solid line 
includes all contributions to T up to next-to-next-to-leading order in the non-relativistic 
expansion. It is obtained from ([8]) assuming that p and p' are parallel to each other. 
The constant dotted blue line gives the leading-order approximation to T. Summing the 
P- or the (momentum-dependent) next-to-next-to-leading order S'-wave contributions to 
the constant S'-wave terms (given by the leading order plus the terms proportional to 5m 
and 5m) yields the dash-dotted red or the dashed black line, respectively. The curves 
are plotted against the relative velocity of the incoming state X1X1 ■ 

(m x o — m x +)/{m x o + ^ x +) or to vertex couplings that vanish in the exact SU{2)l sym- 
metric limit. Similarly, contributions from t-channel exchange amplitudes introduce 
A m factors or coupling factor combinations that lead to vanishing contributions in the 
SU{2) L symmetric theory (case of x Pi waves), or are additionally suppressed (as it is the 
case of 3 Pi and 3 P2 partial-wave configurations) by the masses of t-channel exchanged 
sfermions, since the mass scale of the latter is above 5 TeV in the MSSM scenario consid- 
ered. Consequently, as the initial two particle state in the reaction x\xt ~^ u d consists 
of two wino-like particles with |A. m | ~ 4 • 10 -5 , the 1 P±, 3 P\ and 3 P2 partial waves give 
suppressed contributions to the tree-level annihilation rate. 
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4.3 Example 3: X2X2 ~ ► 

An example of a P-wave dominated process is provided in the left plot of Fig. [2j It 
corresponds to the tree-level X2X2 — ^ h°h° annihilation, wherein S*-wave contributions 
vanish, such that the process is purely P-wave mediated in the non-relativistic regime 
(the coefficient bp is given by 9.94 • 10 _29 cms). The absence of S'-wave contributions 
can be explained by CP and total angular- momentum conservation in the 
reaction^ The CP quantum number of the final two-particle state h°h° is given by 
CP = (— 1) L = (— 1) J , as the total angular momentum of a h°h° state coincides with 
its orbital angular momentum and the parity of such a state is given by P = (— 1) L , 
while its charge conjugation is C = 1. In case of the annihilating xtXa two-particle 
state the J PC quantum numbers are h for a 1 5'o partial-wave configuration and 
for a 3 5'i partial-wave state. Hence, for the xtXa state, CP = —1 is realised in case of 
S- waves for the J = configuration, and CP = +1 for J = 1, which are opposite to 
the CP quantum numbers of a h°h° final state with the same total angular momentum. 
The same reasoning explains the absence of 3 Pi annihilations in any of the processes 
XtXa ~^ X A X B with X A X B = h°h°, h°H°, H H , as the J PC quantum numbers of the 
3 Pi partial-wave configuration of the incoming xtXa states are 1 ++ , hence CP = +1 for 
J = 1. This is opposite to the CP quantum number of the two CP-even Higgs boson 
final state with total angular momentum J — 1. 

Let us finally note that there are also no contributions from 1 Pi partial waves in 
the process shown in the left plot in Fig. [2J This feature is generic to xtXb ~^ X a Xb 
annihilations with identical scalar particles in the final state, XaXb = h°h ,H°H . 
The argument relies on the statistics of the final state identical bosons, and applies 
to all xtXb incoming states and not only to particle-anti-particle states xtXa '■ Bose 
statistics forbids the two identical final state scalars to be in a J = L = 1 state, as the 
corresponding two-particle wave-function for odd total angular momentum J would be 
anti-symmetric. This argument can also be used to explain the absence of the J = 1 3 S'i 
and 3 Pi states in a xtXb ~^ h°h°, H°H° annihilation reaction. 



4.4 Example 4: xTXi W+W -> X2X2 

Let us finally turn to the case of an off-diagonal annihilation rate. The right plot in 
Fig. [2] shows the off-diagonal annihilation rate T associated with the process X1X1 ~^ 
W + W~ — > X2X21 which is relevant, for instance, in the calculation of the Sommerfeld 
enhanced X1X1 ~^ W + W~ and X1X1 ~^ W + W~ (co-) annihilation cross sections. The 
mass splitting between the xt an d xf charginos is given by 324.18 GeV in the MSSM 
scenario considered, which results in rather large mass differences, namely Sm = 8m = 
162.09 GeV. In this case, the Wilson coefficients hi and /12, that are proportional to 

6 The following reasoning applies to all possible xtXa ^ XaXb annihilation reactions with two CP- 
even MSSM Higgs particles in the final state, X A X B = h°h°, h°H ,H°H . Note that CP is conserved 
in these reactions if the mixing matrices in the chargino sector are real, which is the case for the scenario 
we consider. 
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5m and 8m, lead to a 1% positive correction to the constant leading-order rate. This 
positive shift corresponds to the difference between the leading-order approximation to 
the annihilation rate T (first line in flSJ), dotted blue line in the right plot in Fig. [2]), and the 
complete non-relativist ic result for T including next-to-next-to-leading corrections (solid 
blue line) at zero momentum. The corrections induced by the terms proportional to 
5m, 5m turn out to be somewhat smaller than the naive expectation 5m /M = 5m /M = 
2.78%, but represent nevertheless the dominant next-to-next-to-leading order correction 
up to v rc \/c ~ 0.16. For larger relative velocities, the P- and next-to-next-to-leading 
order S-wave terms provide larger contributions to the absorptive part of the XiXi ~^ 
W + W~ — > X2X2 scattering amplitude. This is indicated by the dash-dotted red and 
dashed black curves, which result from the addition of the constant S-wave contributions 
(first two lines in ([8])) and the P-wave contributions (third line in ([8])) or the momentum- 
dependent S'-wave next-to-next-to- leading terms (fourth line in ([8])), respectively. The 
correction to the leading-order T rate, due to the P- and next-to-next-to-leading order 
S'-wave terms amount to a 7% for v rc \/c = 0.4. 

Note that no comparison with public numeric codes providing results for (tree-level) 
XX ~~ X^Xb annihilation rates is available for the off-diagonal annihilation rates. We 
emphasise that the calculation of the partial-wave decomposed off-diagonal annihila- 
tion rates therefore constitutes one of the main results presented in paper I and in this 
work. The relevance of off-diagonal annihilation rates in the calculation of Sommerfeld 
enhanced (co-)annihilation amplitudes in context of the Xi renc abundance calculation 
was in particular pointed out in Sec. 4.2 of I, and will be further investigated in subse- 
quent work [2"T] . 

5 Summary 

With this work we finish the presentation of results associated to the short-distance 
annihilation rates, that are prerequisites for a refined study of Sommerfeld enhance- 
ments in neutralino dark matter (co-)annihilation processes in the MSSM including P- 
and next-to-next-to-leading order S-wave contributions. Our analysis can be applied 
to a set of nearly mass-degenerate non-relativistic neutralino and chargino states with 
masses around the TeV scale and excludes accidental mass-degeneracies with further 
supersymmetric or Higgs particles. 

A factorisation between the short- and long-distance contributions in the pair an- 
nihilation of non-relativistic neutralino and chargino pairs is possible given the large 
separation between the associated scales. Paper I [TH] introduced an effective field the- 
ory set-up (the NRMSSM), that provides the basis for a systematic study of radiative 
corrections in (co-) annihilation processes of non-relativistic neutralinos and charginos, 
applicable to both neutralino DM freeze-out in the Early Universe as well as to neutralino 
pair-annihilations today. In the EFT approach, the tree-level (co-)annihilation rates of 
neutral, single and double charged neutralino/chargino pairs into SM and Higgs two par- 
ticle final states X^Xg, related to the absorptive parts of the 1-loop scattering reaction 
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XeiXe 2 ~ * X^Xb — > Xe 4 Xe 3) are encoded in the absorptive parts of the Wilson coefficients 
of four-fermion operators. As the first step in the construction of the NRMSSM, paper 
I provided the basis of dimension-6 four-fermion operators, which describe leading-order 
S'-wave annihilation processes. The absorptive parts of the corresponding four-fermion 
operators' Wilson coefficients are provided in analytic form in the appendix of paper I. 
In the present work, we extend the results from paper I and provide the operator ba- 
sis for dimension-8 four-fermion operators, contributing at next-to-next-to-leading order 
in the non-relativistic expansion of XeiXe 2 ~>* XaX b — > Xe 4 Xe 3 annihilation rates, and 
present analytic results for the absorptive parts of the corresponding 1 P 1 -wave Wilson 
coefficients as well as the spin-averaged sum of spin-1 P-wave { 3 Pj) Wilson coefficients 
in the appendix |A] An electronic supplement [IE] to this paper contains analytic re- 
sults for all kinematic factors that are needed in the construction of the absorptive part 
of partial-wave separated (next-to-next-to-) leading order S- and P-wave Wilson coeffi- 
cients, relevant for the determination of the 0(v^ cl ) approximation to any (off-) diagonal 
tree-level (co-) annihilation rate. Our results apply to neutralino and chargino states with 
arbitrary composition and include the full mass dependence of the final state SM and 
Higgs particles. 

While there are situations where the main part of the 0(v^ cl ) corrections to a given 
annihilation rate can be attributed to a specific partial wave (for instance when CP- 
conservation or helicity-suppression forbids or suppresses annihilation reactions from 
other partial-wave states), we have presented two examples in Sec. HI where the 0(v^ el ) 
P- and next-to-next-to-leading order S-wave contributions in the annihilation rates are 
roughly of the same order of magnitude and enter with differing signs. A proper treat- 
ment of Sommerfeld enhanced annihilation rates beyond leading order S'-wave annihila- 
tions therefore generally requires the knowledge of each separate partial-wave contribu- 
tion, which is now available with the analytic results given in the appendices of paper I 
and II as well as those collected in the electronic supplement. In particular a numeric ex- 
traction of the 0(v^ el ) contributions in the non-relativistic expansion of the annihilation 
rates without a separation of the different constituting P- and next-to-next-to-leading 
order iS-waves, will, in general, not be sufficient in a rigorous analysis of the Sommerfeld 
effect beyond leading-order S'-waves enhancements. 

In addition, it is worth to stress that our work allows for a consistent treatment of off- 
diagonal annihilation rates, required for the accurate description of Sommerfeld enhanced 
annihilation reactions, as the potential exchange of electroweak gauge-bosons and light 
Higgses prior to the actual annihilation can change the incoming neutralino or chargino 
two-particle state to another nearly on-shell two-particle state. This implies that the 
annihilation process itself is generally described by a non-diagonal hermitian matrix in 
the space of neutralino and chargino two-particle states. Apart from the usual expansion 
in non-relativistic momenta, a consistent treatment of off-diagonal reactions within the 
NRMSSM requires an additional expansion in the mass differences between initial and 
final state particles in the off-diagonal annihilation rates XeiXe 2 ~~ Xj^Xb Xe 4 Xe 3 - 
The respective contributions count as next-to-next-to-leading order in the non-relativistic 
expansion. Consequently, we account for the corresponding set of four-fermion operators 
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in the basis of the dimension-8 four-fermion operators given in this paper and include 
the results for their Wilson coefficients in the electronic supplement. 

With the above results at hand, the study of the long-range effects in the annihilation 
of non-relativistic neutralino and chargino pairs as well as their impact on the neutralino 
relic-abundance calculation in selected examples will be the subject of a forthcoming 
publication |21] . 
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A Absorptive parts of Wilson coefficients of dim- 
ension-8 operators in ££ ann 

We provide in this appendix analytic expressions for the kinematic factors related to the 
absorptive part of the P-wave Wilson coefficients in SC^, Eq. ([3]). The corresponding 
expressions for the kinematic factors of the next-to-next-to-leading order S'-wave Wilson 
coefficients, g{ 2s+1 S s ) and hi( 2s+1 S s ) with s — 0, 1 are quite lengthy. Therefore we have 
collected the latter in an electronic supplement [18] attached to this paper, that also 
contains the kinematic factors associated with the absorptive part of the P-wave Wilson 
coefficients (/( 1 Pi), /( 3 Pj), J = 0,1,2) and those corresponding to the leading-order 
S'-wave Wilson coefficients (/( 1 5 , ), /( 3 Si) in 5Cf^), which were written in the appendix 
of paper I. Details on the nomenclature used in the electronic supplement can be found 
in Appendix IB1 

We aim at the description of Sommerfeld enhanced annihilation rates, and will in a 
forthcoming publication [21] consider the potentials which are responsible for the long- 
range Sommerfeld corrections at leading order. Despite that the leading-order potential 
interactions cannot change the spin of the incoming two-particle state, they depend 
on the spin (s = 0, 1) of the latter. Consequently, as far as our study of Sommerfeld 
enhancements is concerned, the separate knowledge of the different f( 3 Pj) coefficients 
which share the same orbital angular-momentum and spin but different total angular 
momentum ( J = 0, 1, 2), is not needed. It suffices to consider the combination of spin-1 
P-wave Wilson coefficients f( 3 Pj) entering the short-distance part, Eqs. ([8]) and (fTTl) . 

f( 3 Pj) = I f( 3 Po) + f( 3 Px) + ~ /( 3 P 2 ) (15) 
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that will be multiplied by the P-wave Sommerfeld correction factor computed with po- 
tentials for spin-1 scattering states. Hence we give in this appendix analytic expressions 
for the kinematic factors corresponding to /( 1 -Pi) Wilson coefficients as well as the kine- 
matic factors associated with the combination f( 3 Pj). For completeness, the kinematic 
factors for the separate /( 3 Pj), J = 0, 1, 2, Wilson coefficients can be found in the elec- 
tronic attachment [18], together with those of the combination f( 3 Pj). 



A.l Master formula to build the Wilson coefficients 

In paper I we have provided a master formula to obtain the absorptive part f( 2s+1 Lj) of 
a given Wilson coefficient encoding a short-distance annihilation rate of XX states with 
quantum numbers 2s+1 Lj into SM and light Higgs two particle final states XaX b from 
its constituent parts, the kinematic and coupling factors. For the sake of clarity, we write 
the formula here as well and briefly comment on its structure. It reads 

fXe 1 Xe 2 ^X A X B ->-Xe 4 Xe 3 ,2s+l T \ 

J{eie 2 }{e 4 e 3 } \ ^J) 



na 2 I 7 XeiXea-^A^B-^Xea X A X B ,28+1 j \ 



1 



I {^Xe-tXez^XAXB^Xe^Xe-i n (a) X A X B ,2s + l T \ 



a=l n i\,i 2 



a=l n 11,12 / 

Formula ( IT6|) also applies for the next-to-next-to-leading order S'-wave Wilson coefficients 
denoted with g and hi (i = 1,2) in Q, with / being replaced by g or h^. However in 
the discussion that follows, we will generically refer to the absorptive part of any four- 
fermion operator's Wilson coefficient as /. Note that we have used «2 = gf/47r, where 
#2 denotes the SU{2)l gauge coupling. In case of s = 0, the total angular momentum J 
of the 2s+1 Lj state in f|T6l) takes the value J = L, while for s — 1, J — \L — 1\, . . . , L+l. 
The results for f( 2s+1 Lj) at 0(a 2 ) are obtained through matching of the EFT tree-level 
matrix element of four-fermion operators in 5£ ann with the absorptive part of the MSSM 
1-loop XeiXe 2 ~~ Xj±X B -+ Xe 4 Xe 3 scattering amplitude with states XeiXe 2 an d Xe 4 Xe 3 in 
a 2s+1 Lj partial- wave configuration. As tree- level annihilation processes, encoded in the 
absorptive part of the Wilson coefficients at 0(a 2 ), are free from infrared divergences, 
the individual contributions to the Wilson coefficients from exclusive final states X^X^ 
can be given separately. The latter exclusive final state contribution is indicated with 
the label Xe X Xe 2 XaX b -+ Xe 4 Xe 3 o n / in (TTBj) . while the actual absorptive part of the 
Wilson coefficient f( 2s+l Lj) is given by the inclusive annihilation rate, summed over all 
accessible final states. Our results cover separately all possible exclusive SM and light 
Higgs two-particle final states XaX b in neutral (x°X°5 X~X + )i single- charged (x°X + ) 
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X°X~) an d double-charged (x + X + 1 X~X~) chargino and neutralino pair-annihilation re- 
actions, where the XaX b states are conveniently classified to be of vector- vector (VV), 
vector-scalar (VS), scalar-scalar (SS), fermion-fermion (//) or ghost-anti-ghost (riff) 
type, see table 2 in paper I. 

The first line on the right-hand side of ( TT6l) collects all contributions from XeiXe 2 
XaXb —> Xe 4 Xe 3 MSSM selfenergy amplitudes, while the second and third lines give the 
triangle and box amplitudes contributions, respectively Quantities B niii2 , C^ iii2 and 

D"\ li2 in (fl6l) denote the kinematic factors, which encode the 2s+l Lj partial-wave specific 
information on the process. They are obtained from a generic Xe 1 Xe 2 ~~ >* XaX b — > Xe 4 Xe 3 
1-loop scattering reaction with generic external Majorana fermions and generic final 
state particles XaXb of the type VV, VS, SS, ff or r]fj, respectively, and can be 
applied to any XX ~> XaX b XX annihilation reaction with external Majorana or 
Dirac fermions by appropriate construction of the corresponding process-specific coupling 
factors, denoted with lowercase letters (b n i , , c^\ li2 , d„ ili2 ) in f|T6|) . Depending on the 
type of the particles Xa and X B as well as on the topology, there is a fixed number of 
coupling-factor expressions that result from all possible combinations of (axial-)vector or 
(pseudo-) scalar vertex factors at the four vertices of a given MSSM XeiXe 2 ~^ XaX b — > 
Xe 4 Xe 3 1-loop amplitude. These combinations are labelled with the index n in ( [16]) . A 
recipe on how to derive the coupling factors is given in appendix A. 2 of paper I. Note that 
the coupling factors do not depend on the kinematics and hence are the same for Wilson 
coefficients with different quantum numbers 2s+1 Lj. The index a in (TIB]) enumerates the 
expressions related to the four different triangle and box amplitudes (see Figs. 6-8 in 
paper I). Finally, in each of the processes there is a certain set of particle species that 
can be exchanged in the s- or the t-channels of the contributing amplitudes. These are 
labelled with the indices i\ and %2- 



A. 2 Kinematic factors 

Let us first collect from paper I the relevant notation that enters the formulae for the 
kinematic factors. The kinematic factors for a given Xe X Xe 2 XaX b — > Xe 4 Xe 3 scatter- 
ing reaction depend on the external particles' masses, which are rewritten in terms of 
two (in principle) different reference mass scales m, m and two mass differences 5m, 5m 

as 

m ei = m — 5m , m e2 = m — 5m , 

m ei = m + 5m , m e3 = m + 5m , (17) 



such that 



m ei + m e4 m e2 + m 



m = , m 



(■:>, 



2 7 2 1 

6m = m ^ ~ ; 6m = m ^ ~ m e 2 (lg) 
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The mass differences 5m, 5m vanish for the diagonal reactions, XeiXe 2 ~~ Xe 1 Xe 2 ^ while 
they have to be considered as 0(v^ cl ) corrections in the non-relativistic expansion for 
the off-diagonal amplitudes. The convention established by ffTT|) implies that particles 1 
and 4 have masses closer to the reference mass scale m, while particles 2 and 3 share the 
reference scale m. Introducing two distinct mass scales for the particle species allows our 
results for the absorptive part of the Wilson coefficients to cover both the cases of a set 
of particles nearly mass-degenerate with the neutralino LSP (m ~ m) and that of a set 
of non-relativistic hydrogen-like neutralino and chargino systems (m > m or m < m). 
If in a process XiXj ~^ XiXk, the mass mk{m{) is actually closer to the mass m^m.,) 
and the mass scales m and m differ beyond 0(v^ el ), the results for the kinematic factors 
presented below and collected in the electronic attachment cannot directly be used to 
determine the corresponding / expressions, as the mass differences 5m, 5m related to this 
reaction are not necessarily small. However the absorptive part of the Wilson coefficients 
for the reaction with particles 3 and 4 exchanged, XiXj ~* XkXu can be obtained from 
the kinematic factors presented in this work if for that case the particle masses obey (fTT|) 
with mass differences 5m, 5m of 0(v^ cl ). The symmetry relations given in then allow 
to relate the obtained result for the / in XiXj ~^ XkXi rates to the / for the XiXj ~~ ^ XiXk 
reactions. 

We use the hat notation rh a to denote the mass m a rescaled by the mass scale 
M = m + fn, i.e. 



™« = T7 > (19) 



M 



and define the dimensionless quantities 



A m = rh — m 



A AB = m\-m\, 



P = ^l-2(m\ + ml) + A\ B , (20) 

where m A and mg are the masses of the particles Xa and Xb, and for Xa = Xb (3 
corresponds to the relative velocity of the Xa and Xb particle at leading order in the 
expansion in the non-relativistic 3- momenta and the mass differences of the Xe^Q Per- 
forming the same expansion for the single s-channel (gauge or Higgs boson Xj) exchange 
propagators, we obtain the following denominator-structure at leading order: 

P- = • (21) 

The corresponding leading-order expansion of t- and w-channel chargino, neutralino and 
sfermion propagators gives 

PiAB = m m + mf — in rh 2 A — m m B , 



7 For a set of nearly mass-degenerate particles Xe-i > Am will be of the order of the mass differences 
8m and dm, and thus yield terms beyond 0(v^ cl ). The exact dependence on A m is however kept in 
our results, which in particular allows us to cover the case of annihilation reactions in hydrogen-like \X 
systems as well. 
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IB A — J^iAB \A*r*B 



(22) 



The index % in ( 1221) refers to the t-channel exchanged particle species and the labels A 
and B are related to the final state particles Xa and Xb in the actual XX ~^ XaXb 
annihilation reaction. 

It is convenient to write the kinematic factors for the Wilson coefficients of dimension- 
8 operators by pulling out factors of the leading-order propagator and (m m), as well 
as the factor /3 arising from the phase-space integration. For instance, for the kinematic 
factors related to dimension-8 four-fermion operators, that derive from the selfenergy 
topology we define 



tjXaXb /2s+1 r \ _ P r>X A X B (2a+l r 

(m m) z F? Ff 



where the labels %\ and %2 refer to the particle species that are exchanged in the left 
and right s-channel propagator of the selfenergy diagram. As generically either gauge- 
boson (V) or Higgs (S) s-channel exchange occurs in the processes under consideration, 
the combination is given by i\%2 = VV, VS, SV, SS. Note that with respect to the 
definitions in the first paper, there are additional normalisation factors in the prefactor's 
denominator of ( |23i) . Likewise, kinematic factors of dimension-8 Wilson coefficients 
arising from the triangle-topologies are rewritten as 

^(a) X a X b (2s+l t \ P ^(ct)X A X B (2s+l r \ ^10 

U n,hX { L J> ~ Z ^ ~ — G n,uX { L J> Ot—L,Z , 

m m P hAB Fx 

C («)X A X B {2s+lLj) = _P d (a)X A X B (2s+%) a = 3A (24) 

m m P hB A P x 

where the index i\ is now related to the t- or w-channel exchanged particle species, and 
the subscript-index X indicates the type of exchanged particle (X = V, S) in the s- 
channel. Finally, the kinematic factors associated with the box topologies are written in 
the same form as the correspondent expressions related to leading dimension-6 operators: 



n (l)X A X B ( 2s+l r \ _ P f)(l) X A X B /2s+l r 

u n,i x i 2 { L J) - T, B U n,ixn \ L J) 

i%\AB ±i 2 BA 



r ,(2)X A X B( 2s + lr \ P f,(2)X A X B( 2s + lr \ 

L) n,i 1 i 2 { L J) - „ B L> n,i l i 2 { L J) 

-LhAB ±i 2 AB 

r ,(3)X A X B( 2s+l T \ P f,(3)X A X B( 2s+l T \ 

L) n,i 1 i 2 { L J) - ~B B L> n,i l i 2 { L J) 

J BA Jr. i ■ 



BA r i 2 AB 

(4)X A X B( 2s+l T _\ _ P n(4)X A X Bf 2s+l 

\\BA P'i 2 BA 



n (4.)X A X Bf 2s+l T \ _ y n (4j X A X B (2s +\ j \ ( 9 z\ 



In (1251) the indices %\ and % 2 refer to the exchanged particle species in the left and right 
t- and u-channels of the 1-loop box amplitudes, respectively. 
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Finally, let us recall the conventions for the label n established in paper I. Each 
entry for the index n in (12314251) is given by a character string with a length equal to 
the number of vertices that involve fermions in the underlying 1-loop amplitude. In 
case of XaX b = VV, VS, SS or rfq, the string n hence has 2,3 and 4 characters for the 
selfenergy, triangle and box amplitudes, respectively. If X^X B = ff, the string n has 
always 4 characters. The ith element in a string n indicates if the coupling factor at the 
ith vertex of the respective 1-loop amplitude is of vector /scalar (r) or axial/pseudoscalar 
(q) type. The convention to enumerate the vertices in a given diagram is from top to 
bottom and from left to right. Only those kinematic factors with a given label n that 
are non-vanishing are quoted in the following. 



A. 2.1 P-wave kinematic factors for X&X B = VV 

The only non-vanishing kinematic factor ■ in case of l P\ partial-wave reactions is 



given by 

B V qq V vv { l Pi) = ^| (8 /3 2 - 3 A\ B - 27) , (26) 
while for the combined 3 Pj waves the non-vanishing kinematic factors read 

Bl V vv {*Pj) = - ^ (P 2 - 6 A\ B ) , (27) 

B V qq V vv {"Pj) = ± (8 /3 2 - 3 A\ B - 27) , (28) 

B^vsCPj) = Bl v sv CPj) = -\ fh w A m A AB , (29) 

B^ V ss( 3 Pj) = ™w ■ (30) 

In the case XaX b = VV, there are relations among the a = 1(2) and a = 3(4) kinematic 
factors for the triangle and box topologies which are fulfilled for any 2s+1 Lj configura- 



tion (in particular also for the kinematic factors associated with the absorptive part of 
the next-to-next-to-leading order S-w&ve Wilson coefficients, g( 2s+1 S s ) and hi( 2s+1 S s )). 
These can be found in paper I, but are repeated here for completeness: 

M3)VV ( 2s+l j \ _ Ml)VV ( 2s+lr \ I 

M4)VV ( 2s+1t \ _ M2)VV ( 2s+lr \ | 
U n,hV V ^J) — ~ Lj n,i 1 V \ ^J) \A^B , 

^(4)VV/ 2s +lr \ _ M2)VV ( 2s+1t \ I 
j\(Z)VV ( 2s+lj \ _ f)(l)W/2s+lr \ | 
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n(4)VV/2a+lr \ _ p,{2) VV ( 2s+l t \ | foi\ 

The minus sign in the relation for the triangle coefficients C^\^ is a consequence of 
interchanging the two gauge bosons X A and X B at the internal three-gauge boson vertex. 
By virtue of the relations (13T|) . we only need to give the kinematic factors for diagram- 
topologies a = 1,2 for both the cases of triangle and box diagram kinematic factors. 
Starting with the expressions 

for 1 P l partial waves we have 

<5W*) = A» + (A m - 6 Cn„ + 2 A AB ) 

A mm 1^ ^hAB 

'"W (6 A^ A AB - A m (5 (3 2 - 3 A* B ) - 3 A AB ) , (32) 



24 mm 

whereas for the combined 3 Py quantum numbers we find 

C&™v?Pj) = A m A AB + -|!^L. (A m + 2 A^) 

- -4^ (2 A^ A^b - A m (/3 2 - 3 A AB ) + A AB ) , (34) 
8mm 

i(2)VV , 3p x _ /3- 



rO)W /3 n \ 3 gii A _^ 

2 mm ^ 



lAB 

1 



^(5^ 2 -9 + 9Ai-3A, B (A m + A AB )) , (36) 
12 mm 

The coefficients C^]^ S V ( x P\) , corresponding to triangles with a Higgs particle exchanged 
in the s-channel, vanish for all n. The corresponding expressions related to 3 Pj reactions 
read for diagram topologies a = 1, 2 

C%Z('Pj) = ~ - flP^ + ^ (A m A AB + 1) . (38) 

mm t> ^i\AB A mm 

All the remaining non- vanishing kinematic factors C^\^ associated with x Pi and 3 Pj 
scattering reactions with both X = V, S are related to the above expressions by 

Ml)VV /2s+lp \ _ /=y(2) W /2s+1 p \ _ ^(1) VV /2s+l p \ | 

^ggr.iiXV - r -/J — °rg ? ,iiXV ^ J — Wrr.ijJfl ^ J) I WH X ->~ j 

/=y(1)W /2s+lp \ _ /=y(2)VV /2s+lp \ _ A(1)W /2s+l p \i /oq\ 
~~ ^qrq,hX\ r J> ~ ^rqq,ixX\ r J V I ->— ro^ , \,OJ) 
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where these relations hold in particular in case of separate 3 Pj, J = 0, 1, 2, partial-wave 
configurations and hence trivially for the combined 3 Pj waves. Finally, the terms related 
to box diagrams give rise to the following non-vanishing coefficients 

f)Ww ( i p) 



4 (mm) 2 4 (mm) 



^- (A m A AB - 1) 



(A^ (2 - 3 A^) + 3) 



48 (mm) 2 m 12 P% x ab Pi 2 BA 

g 

12 mm P^ab 



(A m A iB + 2m« 2 -l) + {il«5,i 1 Hi 2 }, (40) 



A(2)T/V 
rrrr, i\i 2 \ 1/ 

m, ;i m io m, ; 



4 (mm) 2 4 (mm) 



(A m A AB + 1) 



(A^ (2 /? 2 - 3 A^) - 6 A m A^ - 3) ' 



48 (mm) 2 v m v ^ ' 12 P% x ab P 12 ab 

/3 2 



12 mm -P^ab 

5(1) w f i pl 

rqqr,iii2 V 1/ 



(A m A Ai j - 2 m i2 + 1) + |z! ^ i 2 | , (41) 



4 (mm) 2 48 (mm) 2 

/3 2 



(l2Af n -A^(12-4/3 2 + 3A 2 4B ) + 3) 



(2 A^ + A m (2 m ix - A AB ) - l) 
(2 A m (m^ + m i2 ) + (3 2 + 8 m^m, 

+2 (rh h - fh i2 ) A AB - 2 A 2 AB ) +{ioB,i 1 «! 2 } 1 (42) 



12 mm P^ab 

77-5 5 ^ - * (rra.^ + m i2 j + p~ + 5 m^rn^ 

1^ -r ii^B-r j 2 SA 



n(2)^ flp^ 
rqqr,i\t2 \ 1/ 

rhijhiz /3 s 



4 (mm) 2 12 mm P% x ab 



(2 A^ - A m (2 m n - A AB ) - l) 
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(12 Af n (A m + A AB ) - A 2 m (12 + 4 /3 2 - 3 A 2 ^ 

48 (mmr 



6 A m A AjB + 3) 



02 

+ 19 p p ( 2 A ™ ~ 2 Am ^ + < " 2 A AS ) 



-2 A AB + m i2 - A AB ) - /3 2 + 8 m h m i2 ) 

+ {ii*»i 2 }- (43) 
In case of combined 3 Pj waves we have 



rrrr,i\%2 \ J > 



2 (1-AL) + -^(2AL-A m A^-l) 
2 [mmy 4 [mmy 

+ - (18 A^ + 3 A^ (/3 2 - 2 A 2 ^ - 10) - 4 /3 2 + 6 A 2 B + 12) 

48 [mmy 

2 

+ ~ (5 + A m A AB + 2 (2 m n + m i2 - 2)) 

12 mm P ilAB 

1* -Mi AS -M 2 SA 

+ <-»i 2 } , (44) 



n(2)W / 3p > 



ij^r 2 (l + - Y^kr 2 (2 A 2 ,, + 3 A m A AB + 1) 
2 [mmy 4 [mmy 

+ , *~, 2 (18 A^ + 12 A^ A AB -3A 2 m ((3 2 -2 A 2 AB + 6) + 12 A m A AB 
48 [mmy 

- 4 /3 2 + 6 A 2 AB + 12) 
(3 A^ - A m A AB + 2 (2 m n + m 42 - 2)) 
(9 A 2 m + 18 A m A AB - 12 (/3 2 - 3 m^mj + 9 A 2 AB ) 



VI mm P ilAB 
P 2 



36 Pi lAB Pi 2 AB 
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rqqr,i\i 2 \ 3) 



m i 1 m i 2 a 2 | ^ m h a 

2 [mmy mm 

+ ( 6 + 3 A2 i (/32_ 6 )_8/3 2 + 6Ai B + 6) 

48 [mmy 

B 2 

+ (A^ + A m (2 m h - 4 m l2 - A AB ) - 2) 

12 mm P^ab 

(3 A 2 m -6 A m (m^ + m i2 ) + 2 (/3 2 + 6 rh^rh^) 



12 P% 1 AB Pi 2 BA 



+ 6 (m^-m.J A AB -3 A 2 4B ) + {AdB,!^i 2 } 



rqqr,i\i 2 \ 3) 



1 2 A 2 ^ A r 



2 (mm) 2 '" mm 



(6 A 4 m - 3 A^ (/3 2 + 2) + 12 A m A AB - 8 B 2 + 6 A 2 ^ + 6) 



48 (mm) 2 
12 mm Pi 1 AB 



(3 A 2 t - A m (2 + 4 fh i2 - A A b) - 2) 



B 2 

+ To - !? 5 ( 3 A m ~ 6 A m (m« + - A AB ) — 2 (/3 2 — 6 m^rn^) 

i-4 ±i\AB P 2 AB 



6 (rh h + m i2 ) A AB + 3 A 2 AB ) + jii ^ i 2 } • 



The remaining non-vanishing kinematic factors D^\^ for diagram topologies a 
are related to the expressions given above by 



j\(a)VV /2s+lj \ _ f){a)VV /2s+l r \| 

^qqqqMhy J ) rrrr ,hi 2 \ J I \ m h 2 _5>_m n 2 ' 



(a)VV /2s+lr \ _ f)(a)VV /2s+l r \| 

rrqq,hi 2 \ J) rrrr,iii2\ J }\m,i 2 -^-m i2 1 



qqrr,i\i 2 \ "J) rrrr,iii 2 \ J)\ m i 1 ^— m ii 5 

(a)VV /2s+lr \ _ /2s+lr \i 

qrrq,hi 2 \ ^ J ) ~ ^rqqr^^K ±J J)\m il2 ^>—mi 12 j 
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f\{pt)VV /2s+lr \ _ f)(a)VV ( 2s+l r \| 
^rqrq^ixiiK J ) ~ ^rqqr,^ \ J ) Im^-fr-m^ > 

f)(a)W /2s+lr \ _ n(«)^ /2s+l r \i f 4 o\ 

U qrqr,hi 2 \ ^J) ~ ^rqqrJmK )\m il ^-m. H , {^O) 

where these relations hold for the kinematic factors related to any 2s+l Lj partial-wave 
reaction. 



A. 2. 2 P-wave kinematic factors for X A X B = VS 

The only non-vanishing kinematic factor expression associated with l P\ partial-wave 
reactions and related to selfenergy diagrams reads 

^ 2 

uVS (lp\ _ a 2 

&qq,VV\ ^l) - A m ■ 

In case of combined 3 Pj waves we have 



Bl S vv ePj) = ^ , (50) 



Bl*vs?Pj) = B™sv?Pj) = ^ A m (3 - A AB ) , (51) 



B™ as (?Pj) =^((3 2 -9 + 6A AB - A AB ) . (52) 

The non- vanishing kinematic factors \f related to the four generic triangle topologies 
with gauge-boson exchange V in the single s-channel read 

CiZU 1 ^) = ^S 1 ^ + ^ (Al + A„, - 1 4- A AB ) 4- ^f^- , (53) 
4 mm 8 mm ^ JjiAB 

CSSCPO = <5«™ (^) , (54) 

C^Zer*) = - ^A$, - ^ (AJ„ - A m - 1 + A„) - (55) 
4 mm 8 mm 1^ "ijBi 

CSC-PO = ■ (56) 

In case of combined 3 Pj wave reactions the kinematic factors C^l^f read 

U rrr,hV\ *\7J ~ U rrr,hV\ ^JJ ~ ~ U rqq,hV\ ^) > 

A(3)VS /3p \ _^(4)VS, 3 pN _ ^(3)F5 ,l p v, ,ro\ 

2 mm A mm o 
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M3)VS , 3p ^ 



mw m^ 



2 mm A mm 



^ (Am A^B — 1J — 



6P« 



(60) 



Turning to C^\^ S S factors we find that all kinematic factors corresponding to the 1 Pi 
configuration vanish. Kinematic factors C^\^ S S in combined 3 Pj partial-wave reactions 



read 



C 



(1)VS /3 



^ (3 _ A 



AB 



mm 
1 



(A 



i-3 + 2 m^ t 



16 mm 



(2)VS , 3 p x 
rrr, iiS \ i7/ 



^(3)vs (3p 



/^(1)VS /3p \ 
°7Tr,iiSV "OV ' 

(3 - A AB ) + 

8 mm 
1 



24 Pi 1 AB 

(4 - Aab) A ab + (A m + A m )(3 - A AB )) , (61) 

(62) 



/3 2 



24 R 



i\BA 



(A r 



3 — 2 m,\ 



^ (/3 2 - 3 + (4 - Aab) A ab + (A m - A m )(3 - A AB )) , (63) 
16 mm 



/^(4) VS /3 p \ 



M3) VS ,3 p x 
Wrr.iiSV -OJ 



(64) 



The additional non- vanishing kinematic factor expressions C^Yx whh -X" = V, 5 are 
related to the above given expressions via 



r<X)vs 

^qqr,nX [ 


: 2s+1 p.j) 


_ M2)VS 
— ^rqq,hX [ 




M3)VS 


,2s+l p.j) 


_ M<i) vs 

~ ( - y rqq,i 1 X< 


'2s+lpj 


M1)VS 
^qrq,hX ] 


: 2s+i pj) 


_ ^{2)VS 
— ^qrq^jX 1 


'2s+lpj 


A(3)VS 
^qrq,hX [ 


: 2s+1 pj) 


~ ( - J qrq,i 1 X< 


'2s+lpj 



_ /nr(l) VS ,2s+l n \| 

u rrr,nX\ ^ Im^ -s—m^ > 



(3)VS ,2 8 +l 



rrr, i\X 



P 



J )\rhi 1 -+—mi 1 ' 



C^ S Y ( 2s+1 Pj)\rn, 



J rqq, i\X 

) VS 
'rqq,i\X 



ci 3 l V ; S v ( 2s+1 Pj)U 



(65) 



The above relations hold for kinematic factors C^l^ ( 2s+1 Lj) associated with the ab- 
sorptive part of Wilson coefficients f( 2s+1 Lj) and g( 2s+1 S s ) in and 5Cf^. 
Finally, kinematic factors D^\- for l P\ partial wave reactions are given by 



f)(a)VS /lni 
rrrr,iii2\ 1/ 



/3 5 



24 [mm 



D (i)vs ,i p , 



A 



A 



m, ; 



mi, 



A r 



m, ; 



m, 



[mm] 



,)) 
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~ *^. a (3 Al - 9 + 6A AB ) 
4 [mmy 48 [mmy 

-f3 2 + 3 - 3 (2 - A AB ) A AB ) 

/? 2 



(A m (A m + 2 m l2 - 2) - 1 - 2 
(A m (A m +2m il +2)-l + 2 m i2 ) 
(A m (A m + 2^+2 m ia ) + /3 2 + 4 iri^m, 



24 mm P^ab 
f 

24 mm P 2jBj4 

■t- z -t- z m i2 j -hp -t- 1 m ix Uu l2 
-(2 m n - 2 m i2 + A AB ) A AB ) , (66) 



rqqr,iii2\ 1/ 

= mnA T (A 2 + A m - 1 + A AB ) + mil?? !! 2 A 2 
8 [mmy 8 [mmy 

+ - T^Ti (3(A 2 ra + A m -l)(A 2 ra + A m -l + 2 A Aj? )-/3 2 + 3 A 2 ^) 
96 [mmy 

B 2 

+ (A m (A m + 2 m i2 ) - 2m ir 1) 

24 mm P ilAB 



/3 



2 



24 Pi 1 ABPi 2 AB 



(A m + 2 (m ix + m i2 + A AB )) A m - /3 2 + 4 
+ (2 (m n +m i2 ) + Aab) A ab ) + {*i <+ i 2 } , (67) 



^rqqr,hi 2 ( ^ >1 ) ~ ^rqqr,i!i 2 ( Pl ) I -B- mi 2 > (68) 

fS(4)VS zips ft(2)VS /lpNi , fiq x 

rqqr,i\i 2 \ 1/ rqqr,iii2\ Ulm«m, >~ rfij, 2 ' V Ul ^/ 

The kinematic factor expressions D^] ii2 related to 3 Pj partial-wave reactions read 
D (1)F5 . ( 3 P.) 

A ^ ((^ - ( A m + A ab - 1) - (m n + m t2 ) A m ) - ^gg- A 2 



m 
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+ ^T^ ((3 A^ + 4/3 2 -9 + 6 A AB )A^-3(/3 2 -(l-A AB ) 2 )) 
48 (mmy 

2 

((A m + 2 (2 m n + m i2 + A AB + l))A m + 2 m h + 1) 



24 ram P hAB 
g 2 

24 mm P; 2 £L4 



((A m - 2 (2 m i2 + mi, + A AB + l))A m + 2 m i2 + 1) 



/3 



2 



(3A m (A )) - 12 m^m^ + 



2 



12 -Pi! AS Pi 2 BA 

-3 Aab (2 (m*, + m i2 ) + A Aj r)) 



D (2)T/5 (3 p \ 

rrrr,i\i 2 \ 3) 



^iiA m . . . . m^mjj 



^ (A m (A m + 1) + A^ - 1) - — ^* A^ 

(mm) z 8 (mmy 



96 (mm) 

+ 6 A m (A AB - 1) - 3 (/3 2 - (1 - A AB ) 2 



24 mm P ilAB 



(A m (3 A m + 2 (2 m ix + m i2 + A AB )) - 2 m n - 1) 



+ ox d p f 3 A ™(A m + 2 (mi, + m i2 + A AB )) + 12 m n m l2 - /3 2 

+ 3 A AB (2 (frii! + m i2 ) + A AB )) + jij <b- i 2 } , 



f)(3)V5 , 3 p \ n(l)VS /3p \l 



nW-s (3p N _ nW^ (3p \| 



h(l)VS /3p > 
rqqr,i\i 2 \ 3) 



m i 2 ) A m A A B ~ m h + m i 



4 (mm) A 2 (mm) 



2 



+ ~ T^Ti (A^(/3 2 - 3 A^) + 3) + (A m A AB + 2 m J2 

24 (mmy 12 mm F^ab 



29 



B 2 /3 4 

(A m A AB + 2 m h + 1) + - - , (74) 



12 mm Pi 2 BA 6 P^ab P% 2 ba 



dHYLj'Pj) 



rqqr, i\i 2 ' 



^-(A m A AB + l) + 



2 



4 (mm) 2 4 (mm) 

(A 2 n (p 2 -3A 2 AB )-6A m A AB -3) ' 



48 (mm") 2 v m y ' ~ ; 12 P^P^s 

g 2 

12 mm P ilj4B 



(A m A AB + 2 m i2 + 1) + |i a o z 2 | , (75) 



„(3)VS / 3 p n _ f>(i)vs r 3 P - 

^rqqr^^V r 3> ~ ^ rqqr,hi 2 \ r J ) \m il ^ m i2 i \ ^ ) 

n(4)VS /3p \ . fj(2)VS , 3 p \| (77) 
rqqr,i\i 2 \ 1 J) rqqr,i\i 2 \ 3 J I fh-^m, m il 2 — > —in il 2 " 

Note that relation 072p implies that the denominator structures P- llAB and P% 2BA in the 
kinematic factor corresponding to diagram topology a = 1 have to be replaced by P^ba 
and Pi iAB respectively, in order to arrive at the kinematic factor related to diagram 
topology a = 3. Likewise, in ( I73i) the replacement rule for the kinematic factor for 
diagram-topology a = 2 implies the replacement of Pi lAB and P 2 ab by Pi lBA and Pi 2BA , 
respectively. Similar replacements are needed to obtain the a = 3,4 kinematic factors 
from the a = 1,2 expressions with n = rqqr using (!76l) and (I7T|) . The relations among 
kinematic factors in (j72ti73l) and (j76H77l) also hold for the individual kinematic factors 
related to 3 Pj partial- wave reactions with J = 0, 1, 2. 

The remaining non- vanishing kinematic factors ■ for diagram topologies a = 1, 2 
derive from the above given expressions in the following way: 

fjHVS /2s+lr \ _ (_l\a f)(o<)VS /2s+l r \| 

^ qqqq,i\i 2 \ ^ J ) \ L ) ^ rrrr,hi 2 \ ±J J)\ mi 12 ^f -m il2 j 



j(a)VS /2s+lr \ _ (_T\a+l jj(a)VS /2s+l r \| 
rrgij, iii2\ ^ / \ / rrrr,i\ii\ J )\ 1 ^ l i 2 ~ 



J J) ^ rrrr,iii 2 \ lj J)\ 



f)( a )VS /2s+l j \ _ /_i\q n (a)VS' /2s+l r \| 
-^qrrq^xhV J ) ~ \ L ) ^rqqr^^K Lj J)\ ™i lt2 ~ 



{a)VS ,2s+lj \ _ /_ f)(o) /2s+l r \| 

rqrq,hi 2 \ ^ J I \ 1 ) ^rqqr^^K lj J)\ rn i2 - 



f)( a )VS /2s+l j \ _ _ f)(a)VS (2s+1t \\ ( 7 q\ 
-^qrqr^^K ^ J ) ~ ^rqqr^^K lj J)\ m^-*- -rh il ■ I ' °) 

Similarly, in case of diagram topologies a = 3,4 we find 

n(«)V5 /2s+lr \ _ (i\a f,(a) VS /2s+l r \| 

^ qqqq,iT,i 2 \ ^ J ) \ L ) ^ rrrr,i]_i 2 \ ±J J)\ m il2 ^ -m il2 ; 
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f)(a)VS /2s+lj \ _ f)(a)VS /2s+l r \| 

^rrgg, iii 2 \ / ~ \ 1 ) ^ rrrr,hi 2 \ ±J J)\ rhi 2 -> -m^ ) 

/2s+lr \ _ f)( a ) yS f 2s+1 r M 
1J qqrr,i- l .i 2 \ ~~ ^rrrr, i\i 2 \ lj J)\ m^-* -m il , 

f)(a)VS /2s+lj \ _ f_l\a! fSC")^ 5 *+lr \| 

^qrrq, i\i 2 \ ) \ L ) ^rqqr^inK Lj J)\ m il2 -> ~m il2 i 

f){a)VS /2s+lr \ _ (_i\a f\(a)VS (2s+l j \| 

^ rqrq,i\i 2 \ / v x / ^ rqqr,iii 2 \ ±J J)\ rhi 2 -> -mi 2 ) 

n(a)V5 /2s+lr \ _ f)( a ) yS fWf ^1 ^70^ 

-^(jnjr, iii 2 V J > ~ *^ rqqr,i\ii\ ±J J)\ m^-* -m^ ■ 

The relations in (J75]) - (J79j) are valid for D^]^ ( 2s+l Lj) expressions related to any 2s+1 Lj 
partial-wave reaction. 

A. 2.3 P-wave kinematic factors for XaXb = SS 

In case of XaXb = SS the only non- vanishing kinematic factor B^ s iii2 in 1 P\ partial- 
wave scattering reactions reads 



?SS /l n \ ' A 2 



while the corresponding kinematic factors for combined 3 Pj reactions read 



B^wCPj) = \ Al A AB , 



B S qq S w( ip j) = y > (82) 
BgvsCPj) = B?r S sv( 3 Pj) = % A m A AB , (83) 

B^ssfPj) = • (84) 

The kinematic factors for diagram topologies a = 3(4) and a = 1(2) obey in both the 
cases of triangle and box diagrams certain relations, 

M3)SS f 2s+l T \ _ MX)SS,2s+1t \ I 
U n,iiVl L J) — ~ °n,iiVV L J) \A*+B , 

M4)SS ( 2s+1t \ _ ^(2)5S/2s+lr \ I 

°n,iiVV ^J) — ~ °n,iiVV L J) \A<*B , 

M3)SS ( 2s+lr \ _ Ml)SS ( 2s+l t \\ 
U nMi 2 \ L J) ~ U n,i 1 i 2 \ L J) Uf+B 7 

L> n,i 1 i 2 I L J) - U n ) i l i 2 \ L J> \A^B , (85 J 
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that generically apply for the respective kinematic factors related to a given 2s+1 Lj 
partial- wave configuration, including kinematic factors related to coefficients g( 2s+1 S s ) 
and hi( 2s+1 S s ) (see also Eq. (Ill) in paper I). 

In case of 1 P\ waves we find the following expressions for kinematic factors C^\^y 
and diagram topologies a = 1,2: 

ClZU lp J =C^(^i) = -^A^--|^(A m + 2m il+ A^) . (86) 

24 mm ^ ±%\AB 

In case of combined 3 Pj reactions the corresponding expressions read 

/=*(!) /3 p \ M2)SS , 3 p \ 



rrr, i\ 



= --%A ra A, B -^(A ra A A6 + l) + ^, (87) 
4 mm 8 mm 1^ -t^ab 

Ml)SS ,3p \ _M2)SS / 3p v _ /oon 

12 mm 

Turning to kinematic factors C^^f with a = 1 , 2 we find 



rrr, i 



= ~ ~ z ^ (A m A AB + 1) - — ^ m*! + — — m w , (89) 

8 mm 4 mm 1^ ^ab 

and, as in the case of leading-order 1 S'o and 3 Si kinematic factors (see Eq.(115) in paper 
I), the remaining non- vanishing expressions for C^\ S x with both X = V, S and a = 1, 2 
that are associated with and 3 Pj- (as well as the separate 3 Pj, J = 0, 1, 2) partial- 
wave configurations, derive from the above given expressions in the following way: 

Ml)SS { 2s+lp\ _ M2)SS (2s+lp\ _ _ /y(l)SS ( 2s+l p \ I 

^qqr,hX\ r J) ~ ^ J rqq,i 1 X\ r J) ~ ^rrr,hX\ r J) Im^ -^-m^ , 

Ml)SS { 2s+lp\ _^(2)SS f 2s+lp\ _ _ /nr(l)SS /2s+l p \ | _ / Qn \ 

^grg.iiXl J ~~ ^qr g ,uXl r J) ~ ^rqq,hX\ r J ) \m il ^-m ll ■ 

Finally, the box-diagram related kinematic factors D^\^ for diagram topologies a = 1, 2 
are given by 

A — ( Pl) - " 9^]^ + 24mmP nAB (Am + 2 ^ + Aab) 

(Al + 2A m (m h +m i2 ) 



24 ^ iABp 2 BA 

+(2 + A AB ) (2 m i2 - A AB )) 
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f)( 2 ) ss (1 p i 



/3 2 Am /3 2 A r 



+ 



96 (mm) 2 24 mm P^ab 
P 2 



(A m + 2 m ix + A AB ) 



24 ^ abPi 2 ab 



(A 2 m + 2 A m (m n + fh i2 + A AB ) 



+ (2 m^ + A AB ) (2 m i2 + A AB )) 



(91) 



For the combined 3 P<7 reactions we have 



jjWss (3p \ 

rrrr,i\i2\ 3) 



ni 



[mm) 



\mm) 



(A m A 



+ 



£ 2 



(A m A Aj B - 2 m i2 - 1) 



24 mm P hAB 

(A m A^ B - 1) + \A P ,i a i 2 | 

777 . 777 I ^ I. J 



24 P 14R P 



AB-M2BA 



32 (mm) 



f)(2)SS (3p x 



rrrr, iii2 ' 



(mm) 2 



{mm} 



^- (A m A AB + 1) 



/3 5 



24 mm P^ab 
1 



(A m A 



AB 



2 m i9 + 1) + 



24 fi, 4 Rp 



il 



+ 7^737^7 (A m A AB + l) 2 + {zi^z 2 } , 
32 (mm) 2 

6 2 



Tqqr,i\i2\ 3) 



1 



24 (mm 



(92) 



(93) 



(94) 
(95) 



The remaining non-vanishing kinematic factors can be related to the above given ex- 
pressions by making use of the following relations among D^]f i2 kinematic factors with 
different labels n: 



j\(a)SS /2s+lr \ 
^ qqqq, iii 2 \ J ) 

p,(a)SS /2s+lr x 
rrqq, i\i2 \ Jj 

pXa)SS /2s+lr x 



jj(a)SS /2s + lr 



rrrr, i\%2 1 



n (a)SS /2s+lr \i 
^ rrrr, hi2\ lj J) \ m i2 

{a)SS /2s+l 



-m io j 



rrrr, i\%2 



jj(a)SS /2s+l 
qrrq, i\%2 \ 



Lj ) =Di± S Lj 2s+1 L 



rqqr, i\i2 1 
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f)(a)SS /2s+l r \ _ n(a)SS /2s+l r \| 



rqrq, ii%2\ J / rqqr, 



f)( a )SS (2s+l t \ _ _ f){at)SS /2s+lr \| /q^n 
^ qrqr,%xi%\ n J I ~ ^rqqr^inX ±J J)\ rhi^ -m^ ■ \ JV J 

Note that these relations hold among the kinematic factors associated with any of the 
Wilson coefficients /( 2s+1 Lj), g( 2s+1 S s ) and hi( 2s+1 S s ). 



A. 2.4 P-wave kinematic factors for X A X B = ff 

The relevant kinematic factors Bj i • , related to the selfenergy diagram topology with 
a fermion-fermion final state, read 
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BiLvvQPi) = ^ (/3 2 + 3 - 12 m A fh B - 3 A 2 AB ) , (97) 
for the 1 P\ partial-wave configuration, and 

B f J rr> vvC P j) = -^(/3 2 -l + 4 m A m B + A 2 AB ) , (98) 

B f rr f rr ,vs( 3 Pj) = B f J rr , sv ( 3 Pj) = -^p (m A - m B - (m A + m B ) A AB ) , (99) 

B f rL,ss( 3 Pj) = \ (0 2 + 1 - 4 m A m B - A^) , (100) 

B f q Lvv( 3 Pj) = I (P 2 + 3 - 12 m A m B - 3 A 2 AB ) , (101) 

for the 3 Pj case. In the case that the s-channel exchanged particles are of the same 
type (ii«2 = ^V 7 , SS), the additional non- vanishing kinematic factors are related to the 
expressions (!97I) - fll0ip as 

B'rqgr^i^i Pj) = ^nrr, 1112 ( Pj)\m A m B ^-m A m B j (102) 
Bqrrq^i^i Pj) = Aj<j<j<j, ii «2 ( Pj)\m A m B ^-m A fh B > (103) 

where the notation for the replacement rule applies to the term m A m B , but all other 
occurrences of m A or m B shall be left untouched. Similarly, in case of s-channel particles 
of different type = VS, SV), the additional non- vanishing B^ ili2 terms are given 
by 

Brqqr,i\i-X Pj) ~ ~~ Brrrr,i\i-X Pj) \m A ^-fh A i (104) 
Bqrrq,hi 2 ( + Pj) = ~ Bqqqq,hi 2 ( + Pj)\m A ^-fh A ■ (105) 

There are relations among kinematic factors for diagram topologies a = 3(4) and diagram 
topologies a = 1(2) for both the cases of box and triangle diagrams, that are given by 
(X = V,S) 

r f(3)///2s+lr \ _n(l)///2»+lr \\ 
G n,uxl L J) - U n,hx{ L J) \A<+B , 
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n (3)///2s+lr N 
n (4)///2s+l t \ 



Y(2)//^2s+lr \ I 



D 



(l)///2s+l 



71, «1« 2 1 



'J J U-H-B 



(2) ///2s+l 



L 



J J U-H-B 



(106) 



Note that these relations are valid among kinematic factors associated with any 2s+1 Lj 
partial wave (in particular also for g( 2s+1 S s ) and hi( 2s+1 S s ) associated kinematic factors). 
The structures C^lfy ( 2s+1 Pj) that we obtain for diagram topologies a = 1,2 read 



Ma)}} /ip 1 



A r 



48 mm 



^6 (m A + m B )A AB - 6 (m4 - m s ) 
- A m (f3 2 + 3-12 m A m B - 3 A 2 AB \ 



/3 2 A r 



12 P 



(m A -m B - A AB ) 



i\AB 



and, for the combined 3 Pj partial-wave reactions, 



Ma)}} (3p \ 

^rrrr.iiXK r J) 



16 mm 



(-2 (fh A + m B )A AB + 2 (m A - m B 
+ A m (/3 2 -1 + 4 m A m B + A 2 AB 



P 2 A n 
12 P 



(m A -m B + A AB ) 



i\AB 



fiWff (3p 
rqrq,iiX\ %. 



J) 



(6 (m A + m B - (m A - m B ) A AB ) A„ 

m, m. V 

+ (3 2 + 3 + 12 m A m B - 3 A\ B \ 
(l-m A + rh B ) , 



24 mm 



6 PjAB 

1)// 
^rqrq, i\X 



The following relations for the additional non-vanishing C^jM hold: 



/=.(!) // 
rqrq, iiV 


CPi) 


_/=V(2)// , 
qrqr, iiV* 






l Px) 


mB—t—mB i 


M<*) }} 
qrrq, i\V 


CPi) 


_Ma)ff . 

^qqqq,ilV [ 


a A) 








MVff 

rrqq, i\V 


CPi) 


_M*)ff , 


M) 


. _ / 


'Pl)\ 


m A,B ^ m A,B ■ 


MVff ( 

qqrr, iiV\ 


3 Pj) 


_MVff , 
rrqq, i\V ' 


: 3 ^: 


, . . _ 

1 ^rrrr, iiV 


CPj 


) \mA,B^-mA,B 



(107) 



(108) 



(109) 
(110) 
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ri( a )ff /3n \ _ /3p \| 

Turning to the expressions C^^g, we find that all kinematic factors in case of 1 P\ 
reactions vanish, as it has to be due to total angular-momentum conservation. The 
non-vanishing kinematic factors in combined 3 Pj reactions read (a = 1,2) 



C^JlsfPj) = ^ ( 2 {-rn A + fh B + (m A + m B ) A AB ) A r 

16 mm v - 

+ (3 2 + 1 - 4 m A m B - A 2 AB ^j 



12 Pi 1 AB 

V(l)// /3p \ _n(2)ff /3 n \ /=<(!)// /3 



£ 2 

(l + m A + m B ) , (112) 



^qqrr^jSK r J) ~ Ly rrqq,i 1 S\ r J) ~ ^rrrr,hS\ r J )\rn A , B ->-m AtB ■ \ L10 J 

In case of l P\ partial- wave reactions the kinematic factors D^\i_i{ f° r oi — 1,2 read 
/S(i)// (ip\ 

rrrr, \ 1/ 

((3 2 (1 + A 2 J - (3 - 12 m A m B - 3 A 2 AB ) (1 - A 2 ^ 



384 (mm) 1 

/3 2 / 

+ , ~ [l + m A + m B + A m (m A - m B + Aab, 



48 mm P ilAB 

48 Pi lAB Pi 2 BA 



[l + 2m A + A AB ) (1 + 2 m B — A 



AB 



+ {A*-*- B,h *-H 2 } , (114) 

f)(2)// /ln\ 
rrrr, \ 1/ 

— * (> (A 2 . - 1) + (3 - 12 m A m B - 3 A AB ) (A 2 ., + 1) 
384 (mm) 2 v 

- 12 A m (m A -rh B - (m A + m B ) A AB 

f 

48 mm -P^ab 

48 -F^AB-F^AB 
(lp 

rrqq, V 1 



(l + m A + m B - A m (m A - m B + A AB )\ 
+ aq id td (1 + 2 m A + A AB ) (1+2%- A AB ) + hi <-» i 2 , (115) 
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192 (mm 



\ f/3 2 — 3 + 3 A AB + 12 A m (fh A - rh B - (fh A + m B ) A A b) 

m. m. I V 



A 2 m ((5 2 + 3 + 24 m A m B - 3 A 



AS, 



/3 2 / \ 

+ - — 1 + m A + m B - A m (m A - m B + A AS ) 

48 mm P^ab v ' 



(3 2 ( \ 
+ , ^ ^ 1 - m A - m B - A m (m A - m B - A AB ) 

AH m. m. P: n a \ / 

(/3 2 -4m A m B ) , (116) 



AS mm P i2 BA 

24 Pi 1 ABPi 2 BA 
,(2)ff (I 



rrqq, i\i 2 \ 1/ 

^3- (/3 2 - 3 + 3 A 2 ^ + A^ (/3 2 + 3 + 24 m A m B -3A^ 



192 (mm) 

P 2 



(^1 + fh A + m B + A rn (rh A - m B + A AB ; 



AS mm P h AB 
/3 2 ( 

1 1 - mA - m B - A m (m A -rh B - A 

48 mm P 2 ab ^ 



AS, 



+ o, P ^ B - 4 ^m fl ) • (117) 



24 P% a abP%iAb 



3 i 



The corresponding expressions for combined Pj- partial-wave reactions are 

D^/ri^Pj) = ~ 3 P^i/Pi) + ^ .1' 2 (1 + A 2 J , (118) 

24 (mm) 1 

D?J r f r ? ili2 ( 3 Pj) = 3 D?ll! ili2 { l P x ) + — ^ (1 " K) , (H9) 

24 (mm) z 

fSW// /3 r> \ _ f,(i)// /i p) fh A fh B (3 2 

2 (mm) z o r^ABr^BA 

f)(2)// /3 n \ . f)( 2 )// fip<\ ™ArriB ft 2 ^ ^ , s 

2 (mm) 2 o Pi 1 ABPi 2 AB 

The following relations can be used to obtain the remaining non-vanishing D^\^/ 2 expres- 
sions in case of diagram topology a = 1. Note that they hold for any 2s+1 Lj partial-wave 
configuration. 

f)(l)ff (2s+1t \ _ f)(l)ff (2s+1t \i 

qqqq,ili2\ J > ^rrrr, iii 2 \ J ) \mA,B^-mA,B i 
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L> (1)// . 

qqrr, i\%2 




= D (1)ff - ■ 

rrqq, %\%2 


C s+1 Lj) 




rqqr, i\%2 




= D (1)ff . • 

rrqq, %\%2 


C s+l Lj) 


\m A ^f-m A 1 


D mff . 

qrrq, 1122 




= D Wff . ■ 

rrqq, ii«2 


C s+l Lj) 


|mB— m B j 


rqrq, i±t2 


'2s+l£ \ 


= J D {1)// • 


,2s+l L.j) 


Im^— >—rriA 1 


qrqr, 1122 




= 5 (1)// ■ 


: 2s+l Lj) 


|m s — >-m s • 



In case of diagram topology a = 2 



f)(2)// /2s+lr 
^qqqq^ixiiy J 

f)(2)// /2s+lr 
qqrr,%\%2\ J 

f)(2)// /2s+lr 
rqrq, i\i2 \ ^ 

f)(2)// /2s+lr 
qrqr, iii 2 \ ^ 

5(2)// /2s+lr 

f,(2)// /2s+lr 
qrrq,iii2\ J 



analogous relations exist: 



(2)// /2s+l 



rrrr, ii«2 1 



(2)// , 2s +l 
rrqq, iii2\ 



, f 2s+1 £ 



rrqq, ii«2 x 

(2)// , 2s +l 
rrqq, ii22 \ 

(2)// / 2s +l; 
rrrr, i\i2 \ 



J J IrriA— r«A ) 



J J )\ra A -^t—m A > 



f)(2)// /2s+lr \| 



(122) 



(123) 



A. 2. 5 P-wave kinematic factors for XaX b = rjrj 

The use of Feynman gauge for our computation of the absorptive parts of the Wilson 
coefficients requires to consider unphysical particles in the final states, such as pseudo- 
Goldstone Higgs bosons and ghosts. While the results for final states with pseudo- 
Goldstone Higgses can be obtained from the VS and SS kinematic factors and cor- 
responding coupling structures, the ghosts constitute a different class (rjfj). In order 
to properly construct the coupling factors that go along the kinematic factors B 7 ^ 1 ^ 
presented below, we refer the reader to the rules set up in section A. 3. 5 of paper I. 

For l P\ partial- wave processes, there is only one non- vanishing kinematic factor with 
ghosts in the final state: 



48 m 



J qq,VV 

The corresponding kinematic factors in combined 3 Pj partial-wave processes read 

A 2 



B rr,w( 3]D j) 



16 



24 



m w 



l + A 



AB, 



(124) 

(125) 
(126) 
(127) 
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m w 



A m (1 



m 



(128) 
(129) 



B Notation for the kinematic factors in the elec- 
tronic supplement 

The analytic expressions for the kinematic factors needed to construct the absorptive 
part of the Wilson coefficients up to next-to-next-to-leading order have been stored in 
the Mathematica package attached to this paper. They can be loaded into a Mathematica 
session using the command 



<< kinf actors 



The introductory text in the file explains in detail the notation used for the kinematic 
factors, which we summarise in Tab. |2J 



Kinematic factor 



Name in electronic supplement 



Coupling string n 



nX A X B (2s+l j \ 



Btilde["(ili2)XAXB", n, {2s+l}~L_J] 



il,i2 = V,S 



"rr " , "qq" (X A X B = VV, VS, SS, GG) 



rrrr , rqqr 



"qrrq'V'qqqq" 



(X A X B = //) 



ft{a) X A X B (2s+l t \ 
U n,i 1 X \ L -T) 



Ctilde [alpha, "(X)XAXB", n, {2s+l}"L_J] 



X = V,S 



rrr , qqr 
"rqq'V'qrq 
"rrrr" , "qqqq" , "rrqq 
"qqrr" , "rqqr" , "qrrq 
"rqrq" , "qrqr" 



(X A X B = VV, VS, SS) 



\ (X A X B = //) 



XaXb (2s+1 



( 2s+1 Lj) 



Dtilde [alpha, "XAXB", n, {2s+l}~L_J] 



rrrr , qqqq , rrqq , qqrr 
"rqqr" , "qrrq" , "rqrq" , "qrqr" 



Table 2: Notation for the kinematic factors used in the Mathematica package 

The argument XAXB inside the kinematic factors in Tab. [2] can be given the values 

XAXB = VV , VS , SS , f f , GG 
depending on the type of particles in the final state. The partial-wave configuration 



39 



Quantity 


Name in electronic supplement 




mil, miz 


m A, m B 


m A ml? 

m, mis 




mWr 


A 


um 


A - „ 
t-±AB 




ps ps 


Psil, Psi2 


ps 


PsX 


PixAB-, Pi 2 AB 


Ptil[A,B], Pti2[A,B] 



Table 3: Equivalence between the variables in the kinematic factors introduced in Ap- 
pendix |X] and the corresponding names in the Mathematica package. 



{2s+l}~L_J is specified by one of the following strings: 

"1S0","3S1", for the leading-order S-wave coefficients 

"1P1" , "3P0" , "3P1" , "3P2" , "3PJ" , for the P-wave coefficients 
{2s+l}~L_J = < „ , , 

"lS0,p2","3Sl,p2", for the g( 2s+1 S s ) coefficients 

^ "lS0,dm","lS0,dmbar","3Sl,dm","3Sl ) dmbar", for hi ( 2s+1 S s ) coefficients 

The argument alpha inside the kinematic factors Ctilde and Dtilde in Tab. [2] can get 
as input 

alpha = 1,2,3,4 

referring to our enumeration scheme for the respective four triangle and box topologies. 
Finally, the equivalence between the mass variables and propagator structures introduced 
in Appendix |A] that enter the expressions for the kinematic factors and the corresponding 
names in the Mathematica package are collected in Tab. |3j 
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